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Abstract 


—  The  research  efforts  of  the  research  associates  supported 
by  this  program  are  described.  These  studies  are  divided  into 
two  groups:  I.  Continental  size  arrays  and  their  application, 
and  II.  Theoretical  seismograms  arid  earth  structure. 

A  continental  size  array  can  be  realized  by  considering 
individual  seismic  stations  as  array  elements^  A  section  of 
North  America,  where  LRSM  stations  exist  in  sufficient. density , 
is  considered  as  a  very  large  aperture  array  to  supplement  LASA. 
Travel  times  and  dT/dA  values  for  core  phases  (PKP,  PKIKP)  are 
measured  with  this  array  and  the  results  are  interpreted.  A 
technique  for  direct  measurement  of  dT/dA  is  also  described. 

The  theoretical  seismograms  for  P  and  S  waves  are  computed 
for  a  model  of  point  source  in  layered  medium  and  matched  with 
observed  records!  With  this  technique,  the  structures  are 
determined  for  oceanic  crust  using  marine  refraction  profiles, 
and  for  the  upper  mantle  utilizing  seismograms  from  NTS  nuclear 
explosions . 
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I.  INTRODUCTION 

In  this  report,  the  research  activities  carried  out  during 
the  second  year  of  this  program  are  described.  The  program 
involves  research  associates  and  some  staff  members  from  the 
Department  of  Geology  and  Geophysics  and  the  Lincoln  Laboratory. 
The  main  purpose  of  this  project  is  to  utilize  the  new  tools 
and  techniques  in  seismological  research.  Seismic  arrays  and 
computation  of  theoretical  seismograms  constitute  two  such 
areas  which  are  described  in  this  report. 

The  first  section  of  this  report  is  devoted  to  studies  with 
an  array  of  continental  dimensions.  The  apertures  of  conventional 
arrays  such  ai  LASA  are  too  small  for  seismic  phases  with  very 
high  phase  vexocities.  Extension  of  these  arrays  to  a  thousand 
cr  more  kilometers  is  impractical.  The  technique  that  is 
utilized  here  is  to  treat  a  group  of  seismic  stations  as  an 
array  and  to  use  the  time  of  arrival  at  each  station  to  compute 
the  velocity  (dA/dT)  of  selected  phases  with  a  least  squares 
procedure.  This  is  demonstrated  by  utilizing  the  LRSM  station 
network  as  such  an  array  for  the  study  of  core  phases. 

In  the  second  section  of  the  report,  the  importance  of 
synthetic  seismograms  in  understanding  observed  records  and 
determining  earth  structure  is  illustrated.  At  some  epicentral 
distances  where  a  multitude  of  seismic  rays  emerge  from  a  given 
event  (i.e.  where  travel  time  curves  are  duplicated  or  triplicated) 
seismograms  are  extremely  complicated  and  difficult  to  analyze 
in  terms  of  travel  times  alone.  The  Cagniard-de  Hoop  technique 
is  utilized  to  synthesize  seismograms  which  are  compared  with 
the  observed  records.  The  computed  seismograms  are  very 
sensitive  to  velocity  gradients  in  the  crust  and  mantle  and  to 
the  source  time  function.  Thus  the  method  provides  the  means  for 
studying  both  the  structure  and  the  source  parameter  studies. 
Examples  of  these  are  shown  in  Section  III  of  the  report,  where 
the  method  is  generalized  to  spherical  layer  geometry. 
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CONTINENTAL  ARRAY  AND  THEIR  APPLICATIONS 


II. 1  On  the  Structure  of  the  Earth's  Core 
(E.S.  HusebyA  And  M.N.  ToksOz) 


The  existence  of  a  central  core  in  the  earth  was  first 

suggested  by  Wiechert  in  1897  and  that  of  an  inner  core  by 

Lehmann  in  1936.  Since  then,  there  have  been  a  number  of 

studies  and  different  core  models  proposed.  These  models  vary 

from  each  other  and  they  all  fit  different  aspect  of  the 

data  derived  from  core  phases.  The  variations  arise  primarily 

because  of  the  insensitivity  of  travel  time  curves  to  minor 

structural  features  deep  in  the  core,  and  also  from 

difficulties  of  identifying  branches. 

To  demonstrate  these  problems,  it  may  be  helpful  to 

examine  some  of  the  existing  core  models  proposed  by  Jeffreys 

(1939)  ,  Gutenberg  (1958a)  ,  Bolt  (1964) ,  Adams  and  Randall 

(1964)  ,  and  Ergin  (1967).  These  velocity  models  are  shown  in 

Fig,  1.  Important  parameters  like  travel  time  T,  dT/d.  ,  and 
2  2 

d  T/d.  for  the  above  models  were  calculated  using  classical 

ray  theory  (Bullen,  1961)  .  The  results  of  these  calculations 

are  presented  in  Figures  2-3,  where  the  core  phases  are 

identified  by  the  branches  to  which  they  belong.  From  the 

figures  we  can  deduce  several  important  features.  First,  the 

travel  time  T  is  very  insensitive  to  structural  variations  in 

comparison  with  its  first  and  second  derivatives  (dT/dA  and 
2  2 

d  T/d."  )  .  The  values  of  dT/d.  for  core  phases  ars  small,  and 

errors  in  T  could  result  in  large  errors  in  dT/d.  .  As  a  result, 

2  2 

it  is  almost  impossible  to  calculate  reliable  d  T/dt\  values 

from  observed  travel  time  data.  However,  the  wave  amplitude 

2  2 

A  is  a  function  of  d  T/d  (Thirlaway  and  Carpenter,  1966) : 


A  =  4-, 
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Above  B  is  a  constant,  p,  °  V,  ^  Le  L.  denote  respectively, 
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the  density,  velocity  and  angle  of  incidence  at  the  station  and 

focus,  d§/dT  is  the  time  derivative  of  the  wave  potential,  R 

is  earth  radius  and  H  focal  depth.  The  relative  effects  of 

attenuation  on  different  phases  are  neglected.  In  the  above 

equation  all  quantities  are  either  constants  or  slowly  varying 
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functions  with  the  exception  of  d  T/dA  .  The  expression 
2  2 

dT/dA  '  d  T/d A  will  be  referred  to  as  the  geometrical 
spreading  factor. 

Figure  3  shows  very  clearly  that  amplitude  observations 
are  of  crucial  importance  in  core  analysis.  Large  P  wave 
amplitudes  in  the  distance  interval  A,,'1420-1470  have  frequently 
been  reported  (Shahidi,  1968)  and  different  explanations 
have  been  given  for  this  feature.  Basically,  the  above 
phenomenon  is  attributed  to  a  large  geometrical  spreading 
factor,  which  may  be  combined  with  diffraction  effects  near  the 
caustic  B  (Figure  3) .  The  latter  hypothesis  (Jeffreys,  1939) 
favors  the  interpretation  that  the  precursors,  i.e.  waves 
arriving  prior  to  the  P(DF)  phases,  are  real  core  phases  and 
not  diffraction  effects.  Unfortunately,  it  is  very  difficult 
to  verify  diffraction  phenomena  associated  with  a  caustic  in 
this  case,  as  the  local  velocity  variation  is  not  well  known 
and  might  strongly  influence  the  diffraction  pattern.  Another 
difficulty  here  is  the  problem  of  proper  phase  identification 
in  the  140-150*  distance  interval. 

In  recent  years,  much  attention  has  been  given  to  the 
precursor  branches  mentioned  above.  Although  these  waves  are 
relatively  easily  detectable  on  modern  instruments,  the  number 
of  phases,  at  most  3  so  far,  and  their  arrival  times  are 
somewhat  controversial.  An  important  question  here  is  how 
many  core  phase  branches  can  be  observed.  This  question  has 
to  be  raised  as  dT/dA  for  core  phases  vary  continuously  between 
say,  0.0-4.60  sec/deg  and  the  well  established  DF  and  AB 
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branches  between  A  ®143°and  180*  account  for  roughly  two- 
thirds  of  this  variation.  Also,  the  precursor  phases  are 
always  very  weak.  This  is  not  compatible  with  the  values 
of  their  geometrical  spreading  factors.  As  a  result  one 
wonders  whether  there  may  be  another,  quite  different 
explanation  of  these  waves  than  those  advocated  up  to  now. 

Another  property  of  existing  core  models  is  the  small 
2  2 

d  T/dA  values  and  large  dT/dA  values  allocated  for  most  part 
of  the  AB  branch.  The  first  term  gives  a  small  geometrical 
spreading  factor  while  the  latter  requires  a  large  angle  of 
incidence  at  the  mantle-core  boundary.  That  is,  at  large 
distances  the  amplitude  Of  the  P(AB)  would  be  weak  compared  to 
the  P (DF)  phase. 

In  this  study,  in  addition  to  travel  times  we  measured 
dT/dA  and  the  relative  amplitudes  of  the  core  phases  for 
determining  a  structure  which  is  constrained  by  all  these 
quantities. 

DATA  AND  METHODS  OF  ANALYSIS 

To  determine  travel  times  of  core  phases  a  large  number 
of  seismograms  from  North  American  WWSS  and  LRSM  stations  were 
analyzed  (Table  I)  for  fifteen  earthquakes  listed  in  Table  II. 
The  procedure  used  was  to  identify  the  wave  associated  with  a 
certain  branch  and  to  read  the  arrival  time  at  each  station. 

To  determine  dT/dA  for  core  phases  an  array  with  the  proper 
dimensions  is  needed  for  the  measurement  of  such  high  phase 
velocities.  Existing  arrays  such  as  LASA  are  too  small  for 
measuring  phase  velocities  larger  than  25  km/sec. 

We  consider  the  North  American  station  network  as  a  large 
array  and  express  the  travel  time  variation  of  this  array  as 
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where  ^  and  A  are  station  latitude  and  longitude.  Expansion 
of  in  a  Taylor  series  gives 

j  A  ^  i  A*)-*  4-  ^  cU  '  Ao^) ( 3 > 


Neglecting  terms  of  higher  order  than  two,  we  can  write  for  T, 
at  the  i-th  station 

—  \  ^ 

\i  -  a.  +<*.,<}; +  cw}i  +  <M<Ma**s  (4) 

Now  T  is  expressed  as  a  second  order  polynomial  in  §  and  X  . 
From  a  set  of  observations  by  the  method  of  least  squares 
the  coefficients  can  be  determined. 

The  ray  parameter  p  (i.e.  the  gradient  of  the  scalar  T) 
is  obtained  by  differentiating  eq.(4). 

where  a  bar  denotes  vectors ,  The  azimuth  of  the  normal  to  the 


wave-front  is 


A-*;  y  e ^ 


By  taking  the  derivative  of  the  (5)  in  the  direction  of  ray, 
we  obtain  d^T/dfc^ 
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It  should  be  noted  that  this  procedure  is  independent  of 
focal  parameters.  In  fact,  approximate  epicenter  location  is 
easily  obtained  from  the  calculated  azimuth  values  and  a  table 
of  p-  values  as  function  of  distance.  As  stated  before, 
accurate  travel  times  are  of  critical  importance  when  dealing 
with  core  phases.  Therefore,  wherever  possible,  the  records 
were  correlated  with  each  other  to  minimize  reading  errors 
(Jansson  and  Husebye,  1968) ,  also  available  station  corrections, 
namely  those  of  (Chinnery  and  Toksflz,  1967)  and (Cleary  and 
Hales,  1966)  were  incorporated  in  the  analysis.  The  former 
gave  rise  to  the  smallest  standard  error  in  T  and  were  adopted 
in  this  study.  To  reduce  errors  of  this  type,  the  arrays 
formed,  consisted  of  stations  in  the  NE  or  SW  part  of  North 
America.  During  the  dT/d A  calculation,  array  size  and 
distance  range  were  varLed  when  sufficient  data  were  available. 
This  procedure  proved  useful  in  locating  errors.  Since  the 
accuracy  of  dT/dA  decreases  with  distance  from  the  array 
center,  only  the  central  dT/dA  value  was  retained  for  further 
analysis.  The  only  exception  was  the  BC-branch  where  few 
records  were  available.  Accordingly,  despite  the  large  number 
of  records  used,  the  dT/dA  data  obtained  were  few. 

In  the  case  of  core  phases,  the  main  problem  is  still  a 
proper  phase  identification.  The  amplitude  of  these  phases 
play  an  important  role  in  the  identification.  Wherever  possible 
amplitudes  were  measured  and  computed.  To  avoid  the  problems 
of  normalizing  magnitudes  and  varying  crustal  effects  at 
different  stations,  amplitude  ratios  such  as  Amp (AB) /Amp(DF) , 
between  branches  were  used  for  each  event. 
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EXPERIMENTAL  RESULTS 

The  T,  dT/dA  and  amplitude  values  obtained  for  the  main 
branches  of  the  core  phases  are  presented  in  Figures  5-8.  For 
some  distance  interval  our  dl’/dA  values  deviate  significantly 
from  those  previously  reported.  A  more  detailed  presentation 
of  the  data  for  each  branch  is  given  below. 

The  DF-Branch 

The  main  problem  here  is  to  locate  the  branch  endpoint  D. 

The  core  phases  have  been  observed  down  to  about  105°. 

Gutenberg  (1958a)  and  Ergin  (1967?  put  D  around  A  *  125° ,  while 

Jeffreys  (1939)  places  it  at  A  ■  110°.  We  terminate  the  DF- 

branch  at  A =123°,  as  the  amplitude  curve  has  a  pronounced 

maximum  here  (Figure  9) .  If  the  inner  core  boundary  is 

characterized  by  a  sharp  velocity  increase,  then  the  amplitudes 

of  P (DF)  waves  near  the  endpoint  D  should  be  small  due  to  nearly 

2  2 

grazing  incidence  and  small  d  T/dA  values.  In  our  opinion, 
the  waves  observed  between  105*  -  123*  do  not  belong  to  the 
DF-branch  but  are  reflections  from  an  inner  core  boundary. 

Some  records  of  P-waves  in  this  distance  range  are  displayed 
in  Figure  10. 

The  DC-Branch 

This  branch,  sometimes  denoted  as  the  receding  branch  of 
DF,  ranges  from  A*  123*  to  C  at  about  160*  where  dT/dA(V2.45  sec/deg. 
From  what  has  been  said  above,  it  seems  appropriate  to  prolong 
this  branch  down  to  105*,  naming  the  endpoint  D*.  The  obser¬ 
vational  evidence  of  the  existence  of  the  DC-bra.ich  is  meager, 
and  has  been  based  mainly  on  the  requirement  of  t  monotonic 
increase  of  dT/d.'.  between  0.0  and  4.60  for  core  phases.  Subiza 
and  BAth  (1964)  attributed  observations  of  large  wave  amplitudes 
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at  A =130*  -  138*  to  this  branch.  Some  of  our  records 
strongly  support  the  above  observations  as  demonstrated  in 
Figure  10.  Identification  errors  are  possible  since  the  travel 
time  difference  T (DF)  -  T(DC)  would  be  about  2-4  seconds. 

It  should  be  noted  that  only  one  of  the  events  analyzed  gave 
a  dT/d.  value  which  is  expected  for  the  DC  branch.  On  the  other 
hand,  the  assumed  DF -waves  sometimes  give  dT/d A  values 
expected  for  the  DC-branch.  Thus,  the  given  dT/dA  values  for 
the  DC-branch  are  based  on  the  endpoint  values  of  the  DF  and 
BC  branches  and  amplitude  data.  Amplitudes  along  the  DF- 
branch  should  decrease  steadily  toward  the  endpoint  D  due  to 
decreasing  geometrical  spreading  factor  and  energy  transmission 
through  an  inner  core  boundary.  For  large  amplitudes  on  parts 
of  DC,  the  range  of  variation  of  dT/d A  must  be  sufficiently 
large.  Our  value  is  about  0.5  sec/deg  while  corresponding 
values  for  Adams  and  Randall  and  Bolt's  models  are  0.20  and 
0.10  sec/deg  respectively.  Also  their  distance  range  is  about 
15°  larger  than  ours. 

The  Precursor  Branches 

An  interesting  but  confusing  feature  of  core  phases  is 
the  existence  of  some  precursors  before  the  DF-branch.  These 
are  named  GH  and  IJ  branches  (Figure  2)  by  Adams  and  Randall 
(1964)  .  These  forerunners  were  first  attributed  to  diffraction 
effects  near  a  caustic  at  143*  (Bullen  and  Burke-Gaf fney ,  1958) . 
Later  studies,  (Bolt,  1964;  Adams  and  Randall,  1964)  introduced 
a  layered  zone  outside  the  inner  core  as  a  more  satisfactory 
explanation  of  these  waves.  The  precursor  branches  will  not 
be  included  in  our  core  model  for  several  reasons.  First, 
precursors  are  observed  from  about  A  =  14 4 c  and  down  to  about 
105*  and  they  do  not  resemble  ordinary  P-waves  (see  Figures  10 
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and  11) .  They  are  in  general  long  and  oscillatory  wave  trains 
similar  to  guided  waves.  It  is  noteworthy  that  waves  identified 
as  forerunners  on  some  long-period  records  (  A»105#  -  112*) 
give  evidence  of  negative  dispersion.  Second,  the  calculated 
dT/dA  values  partly  overlap  with  those  of  our  BC-branch.  The 
precursor  branches  are  characterized  Ly  very  weak  onsets. 

This  may  explain  the  differences  in  arrival  times  between  Adams 
and  Randall's  IJ-branch  and  Bolt's  FG  branch  which  is  about 
7  sec  at  A-  135°.  Our  own  records  sometimes  give  arrival 
times  relative  to  T{DF)  in  agreement  with  Adams  and  Randall's 
data,  other  times  with  those  of  Bolt,  and  in  a  few  records, 
arrivals  precede  those  of  Adams  and  Randall  by  a  few  seconds. 

The  possibility  of  three  precursor  branches  as  proposed  by 
Ergin  (1967)  cannot  be  excluded.  The  above  points  are 
important,  as  the  precursors  arrive  between  the  two  main  phases 
(AB  and  DF)  at  distances  beyond  140°.  At  these  distances  it 
is  not  possible  to  vary  the  phase  arrival  times  much.  This 
explains  the  large  differences  in  reported  dT/d&  values  for  the 
precursor  branches.  From  the  available  evidence  it  is  not  clear 
that  the  forerunners  are  core  phases  of  ordinary  type  as 
interpreted  by  Bolt  (1964) ,  Adams  and  Randall  (1964) and  Ergin 
(1967) .  If  we  interpret  the  above  precursors  at  least  partly 
as  some  sort  of  guided  waves  (to  be  discussed  later) ,  then 
our  dT/dA  values  will  correspond  to  the  slowness  of  energy 
bursts  in  these  waves.  Also,  such  an  explanation  will  remove 
the  inconsistency  in  the  reported  arrival  times  and  number  of 
phases  identified.  Precursors  are  very  weak  between  A  »114° 
and  A  s  124°,  and  this  may  indicate  that  more  than  one 
excitation  mechanism  exists. 

The  BC-Branch 


The  BC-branch  displayed  in  our  Figure  6  incorporates 
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features  from  different  models,  but  it  does  not  agree  wholly 

with  any  single  model  proposed  previously.  The  branch  we 

have  shown  is  similar  to  that  of  Adams  and  Randall  (1964)  in 

the  sense  that  it  excludes  the  possibility  of  strong  diffraction 

effects  near  B.  Also,  like  Jeffreys  and  Gutenberg's  models, 

dT/d&  values  range  from  about  2.45  -  3.45  sec/deg.  The 

hypothesis  of  strong  diffraction  effects  was  omitted  for 

several  reasons.  First,  observed  dT/dA  values  do  not  favor 

such  features,  the  amplitude  observations  of  P(AB)  waves  require 

2  2 

that  dT(AB)/dA  is  almost  constant  near  B.  The  d  T(BC)/dA 
parameter  is  large  enough  to  explain  the  observed  amplitudes 
along  the  BC-branch.  Theoretical  amplitude  curve  given  in 
Figure  8  is  the  geometrical  spreading  factor  of  the  BC-branch 
relative  to  the  DF~branch.  If  B  were  associated  with  diffraction, 
then  we  should  expect  irregularities  in  signal  shape  and 
amplitude  variation.  However,  no  evidence  of  such  behavior 
of  P (BC)  could  be  detected.  This  check  was  performed  partly 
by  visual  examination  and  partly  by  applying  a  signal-matching 
method.  The  latter  method  also  was  used  unsuccessfully  in 
search  of  the  fourth  P  phase  of  Adams  and  Randall's  and  Bolt's 
core  models,  i.e.  P(IJ)  or  P(GH).  In  fact,  even  the 
identification  of  the  AB-branch  was  questionable  fc  l  Ax  141*-152°  . 
An  exact  location  of  the  branch  endpoints  B  and  C  is  impossible 
due  to  decreasing  amplitudes  toward  these  points.  The  dT/dA 
estimation  was  difficult  as  the  data  available  was  concentrated 
at  Ax  144. -  147°  and  A  «■  150°  -  152°.  Between  A  *141°  and 
144.5°  the  time  lags  between  T (DF) ,  T(DC),  T(BC)  and  T (AB)  are 
small,  preventing  reliable  phase  identification. 

The  AB-Branch 


The  amplitude  observations  in  Figure  8  require  a  slow 
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dT(AB)/dA  increase  all  the  way  from  B  toward  A,  contrary  to 
that  in  other  core  models.  Our  amplitude  data  and  those  of 
Shurbet  (1967)  and  as  well  the  few  dT(AB)/dA  data  available 
confirm  this  hypothesis.  In  the  distance  interval  155°-180#, 
all  T (AB)  and  T (DF)  were  obtained  from  records  of  a  single 
event. 

Travel  time  analysis  for  stations  near  the  antipodal 
point  are  difficult  as  epicentral  distance  may  be  nearly 
constant  while  station  separations  and  azimuths  very  considerably. 
In  other  words,  even  very  small  lateral  velocity  gradients  in 
the  mantle  and  outer  core  could  cause  large  errors  in  the  dT/dA 
values.  Observed  and  calculated  travel  time  data  of  P(DF), 

P(BC)  and  P(AB)  are  shown  in  Figure  5.  A  clear  secondary 
onset  occuring  at  both  the  AB  and  DF  branches  is  attributed  to 
crustal  reflection  or  reverberation  effects.  The  location  of 
the  endpoint  A  and  also  the  corresponding  dT/dA  value  are 
uncertain.  From  our  data  we  cannot  determine  whether  AB-branch 
extends  beyond  180°,  while  Ergin  (1967)  reports  weak  P(AB) 
onsets  out  to  200*. 

A  NEW  CORE  MODEL 

The  dT/d A  data  we  presented  in  the  previous  section  differ 
from  those  of  other  core  models  at  some  distance  intervals. 
Although  there  are  regions  where  observations  are  scarce,  the 
new  results  warrant  interpretation  in  terms  of  a  velocity 
structure  inside  the  core.  The  procedure  we  use  is  a 
simplification  of  that  outlined  by  Bullen  (1961)  which  includes 
stripping  off  of  the  mantle.  Although  the  above  method  is 
straight  forward,  we  sometimes  have  to  resort  to  th'-*  reverse 
procedure,  that  is  fitting  the  preliminary  model  to  the  observed 
T,  dT/d&  and  amplitudes.  The  reason  for  this  is  the  insensitivity 
of  the  above  method  to  small  travel  time  gradients  near  branch 
endpoints . 
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Before  we  proceed  with  a  detailed  discussion  of  the  new 
core  model,  we  should  discuss  briefly  the  accuracy  of  the  data. 
The  standard  deviation  of  the  dT/dA  estimates  are  in  general, 
small,  around  -0.08  sec/deg.  The  records  from  the  LRSM  stations 
are  of  an  excellent  quality,  permitting  reading  of  arrival 
times  to  the  nearest  0.1  sec.  However,  manipulations  with 
array  size  and  station  corrections  indicate  existence  of  biased 
errors  which  may  amount  to  -0.10  sec/deg.  This  figure  is 
likely  to  be  larger  for  the  BC  and  AB  branches.  We  searched 
for  large  errors  in  the  above  dT/dA  values  by  checking  observed 
amplitude  ratios  with  those  computed  on  the  basis  of  our  model. 
The  effect  of  errors  in  dT/dA  on  the  P-velocity  depth 
distribution  is  most  easily  anticipated  at  branch  endpoints. 

For  example,  a  dT/dA  interval  of  0.02  sec/deg  is  equivalent  to 
a  depth  interval  of  about  11  km.  The  corresponding  velocity 
increment  is  small,  about  0.006  km/sec  for  parts  of  the  AB- 
branch.  Just  near  boundaries  where  velocity  discontinuities 
exist,  distoit^on  of  these  boundaries  will  account  for  larger 
velocity  errors.  This  is  demonstrated  in  Figure  1,  which 
shows  that  the  given  core  models  with  exceptions  of  Ergin's 
and  a  small  part  of  Jeffreys"  are  very  similar  despite  large 
differences  in  the  original  data  from  which  they  were  derived. 
Although  the  absolute  values  of  travel  times  are  influenced 
by  focal  parameters  and  the  velocity  distribution  in  the  mantle, 
the  time  differences  between  the  BF,  BC  and  AB  branches  should 
be  insensitive  to  parameters  other  than  the  epicentral  distance. 

We  consider  the  amplitudes  of  core  phases  as  a  strong 
constraint  in  the  interpretation  of  travel  time  and  the 
dT/dA  data.  The  joint  interpretation  of  travel  time,  dT/dA 
and  the  amplitudes  remove  much  of  the  ambiguity  otherwise 
involved.  A  detailed  discussion  of  this  important  point  is 
given  by  Asbel,  et  al.  (1966) .  As  amplitude  observations  of 
short  period  P-waves  exhibit  considerable  scattering  we  only 
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seek  a  broad  agreement  between  the  observed  and  the  calculated 
amplitudes.  To  minimize  source  magnitude  and  crustal  response 
effects  we  use  the  ratios  of  the  amplitudes  of  different 
phases  such  as  A(AB)/A(DF). 

A  factor  of  crucial  importance  in  deducing  a  core  model 

is  the  value  of  dT/dA  at  the  core-mantle  boundary.  For 

example,  commonly  accepted  P-wave  velocity  values  for  the 

lowest  200  km  of  the  mantle  is  around  13.60  km/sec,  which  gives  a 

maximum  dT{AB)/dA  value  of  about  4.45  sec/deg.  However, 

this  is  incompatible  with  our  observations.  One  way  to  avoid 

this  difficulty  is  to  increase  the  dT(AB)/dA  factor  by 

lowering  the  P-velocity  in  the  lowest  part  of  the  mantle.  Our 

final  values  are  shown  in  Figure  14.  Such  a  velocity 

distribution  and  dT/dA  values  of  4.55-4.60  sec/deg  are  in 

general  agreement  with  the  recent  array  measurements  and  the 

resultant  structure  of  the  deepest  part  of  the  mantle  (Johnson, 

et  al. 

1967;  Sachs,  1967;  Fairborn,  /y!96 7)  .  At  the  endpoint  A  of  this 
branch,  the  depth  of  penetration  is  about  2530  km  from  the 
earth's  center.  This  means  that  the  velocity  distribution 
from  this  depth  to  the  core-mantle  boundary  is  obtained  on  the 
basis  of  travel  times.  In  other  words,  with  only  PKP  and 
PKIKP  data  jt  is  quite  permissible  to  vary  the  velocity 
distribution  in  the  outer  core  as  long  as  total  travel  time  is 
satisfied.  We  had  to  assume  a  slightly  steeper  velocity 
gradient  in  the  outer  core  than  that  of  Jeffreys'  model.  This 
was  necessary  to  ensure  reasonable  travel  times  for  P(AB).  The 
dT(AB)/dA  and  amplitude  data  require  an  almost  constant  velocity 
gradient  between  R  =  2500-2150  km. 

As  for  the  BC-branch,  the  effect  of  small  errors  in  the 
velocity  distribution  in  the  AB-depth  interval  will  diminish 
toward  C.  Assuming  that  A(BC)/A(DF)  is  mainly  governed  by  the 
geometrical  spreading  factor,  our  model  gives  a  reasonable  fit 
of  expected  amplitude  ratios  to  those  actually  observed  (and  to 
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those  of  Shurbet,  1967) .  The  calculated  P-wave  velocity 
distribution  for  the  BC  interval  may  be  interpreted  as  a 
transition  from  homogeneous  material  at  B  to  inhomogeneous 
material  at  C.  The  branch  point  B  is  characterized  by  a 
turning  point  in  the  curvature  on  the  velocity  curve.  We 
would  like  to  point  out  that  small  variations  of  the  velocity 
structure  here  would  permit  an  extension  of  the  AB  and  BC- 
branches  toward  shorter  distances. 

As  pointed  out  previously,  hardly  any  reliable  observations 
of  the  P(DC)  phases  are  available.  This  means  that  the  structure 
near  the  inner  core  boundary  may  be  more  complex  than  hereto 
assumed.  A  rapid,  continuous  velocity  increase  in  this  region, 
as  favored  by  Gutenberg,  produces  structures  which  would  not 
permit  core  wave  observations  at  distances  shorter  than  A  0 i  125°. 
The  branch  endpoint  D  is  the  point  where  the  reflected  P  wave 
from  the  inner  core  boundary  obtains  its  critical  value  and  is 
located  somewhat  arbitrarily,  at  A  *  1?3°.  This  distance  was 
chosen  because  at  A  ■  120°  -  124°  the  wave  amplitudes  are 
very  large  and  some  signal  decorrelations  occur.  The  latter 
effect  may  partly  be  due  to  the  dominance  of  high  frequency 
energy  in  the  signals.  Actually,  we  should  expect  some 
pecularities  in  the  P-wave  near  the  critical  angle  of  incidence. 
Another  factor  favoring  location  of  D  at  A  «  123°  is  that  this 
would  allow  a  larger  energy  concentration  in  the  region  A" 110°  - 
123#  as  required  by  amplitude  observations. 

The  DC-branch  is  constituted  of  waves  reflected  from  the 
inner  core  boundary,  and  more  formally  should  be  denoted  as 
PKiKP-waves.  We  find  it  convenient  to  locate  D'  at  &  =  105° 
as  it  is  very  difficult  to  observe  this  wave  for  shorter  distances 
due  to  decreasing  reflection  coefficients  and  interference 
from  other  waves. 

In  the  previous  section  we  discussed  the  reasons  for  our 
preference  to  interpret  the  precursor  branches  as  guided  waves 
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rather  than  core  waves  of  ordinary  type.  As  guided  waves 
travel  time  puts  a  strong  constraint  on  possible  sources  of 
these  waves.  First,  they  cannot  be  associated  only  with  an 
inner  core  boundary  or  the  transition  zone  as  their  travel  times 
do  not  fit  such  paths.  The  outer  core  is  not  a  likely  region. 

The  ray  deflection  at  the  mantle-core  boundary,  would  make 
the  angle  of  incidence  in  this  region  too  small  to  trap  sufficient 
wave  energy.  Focal  mechanism  peculiarities  have  to  be  excluded 
since  the  precursor  observations  are  limited  to  certain 
distance  intervals.  One  possible  explanation  is  propagation  in 
the  mantle-core  transition  zone.  So  far,  only  the  propagation 
of  high  velocity,  low  frequency  diffracted  waves  has  been 
investigated  in  this  zone  (Phinney  and  Alexander,  1966). 

However,  slow,  high  frequency  diffracted  waves  as  reported  here 
may  be  possible,  and  this  problem  is  being  studied  at  this  time. 
The  above  explanation  may  not  be  sufficient  to  account  for  all 
the  features  of  the  precursors.  Similar  excitation  might  exist 
at  the  inner  core  boundary  causing  the  observed  increase  in 
precursor  amplitude  at  distances  beyond  128°.  Some  contribution 
also  may  come  from  ordinary  P-waves  as  a  result  of  possible 
prolongation  of  the  P(AB)  and  P(BC)  branches  due  to  velocity 
peculiarities  at  R~2000  km. 

Our  new  model  for  the  core  and  the  lower  most  part  of  the 
mantle,  is  presented  in  Figure  14.  We  have  introduced  a 
velocity  decrease  in  the  mantle  and  a  gradient  increase  in  the 
outer  part  of  the  core.  These  changes  are  based  entirely  on 
our  analysis  of  PKP  waves  traveling  in  the  lower  parts  of  the 
core.  However,  the  essential  feature  of  the  new  core  model 
is  an  inhomogeneous  zone  extending  out  to  R=  2530  km.  The 
main  evidence  for  the  inhomogeneity  is  the  new  values  of 
dT(AB)/dA.  Although  seemingly  modest,  they  are  responsible 
for  this  extension.  Furthermore,  this  effect  is  not  removable 
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by  small  changes  in  the  velocity  distribution  in  the  outer  part 
of  the  fluid  core.  At  ^.^  2000  km,  corresponding  to  the  branch 
point  B,  the  velocity  gradients  increase  rapidly  over  a  small 
depth  interval.  We  cannot  exclude  the  possibility  of  a  small 
velocity  discontinuity  in  this  region.  For  the  remaining  part 
of  the  core  (R<1900  km),  the  velocity  gradient  is  small, 
positive  or  negative,  which  is  in  good  agreement  with  other 
models.  The  much  discussed  precursor  branches  would  have  a 
modest  influence  on  the  velocity  structure.  Our  calculations 
give  that  the  inner  core  radius  is  around  1250  km,  and  this 
boundary  may  be  preceded  by  some  sort  of  a  transition  zone. 

The  exact  location  of  the  branch  point  D  does  not  strongly 
effect  the  above  structure. 

Assumptions  introduced  when  deriving  the  structure  of  the 
outer  core,  naturally  would  affect  the  calculations  of  the 
inner  core.  In  other  words,  the  core  model  above  might  need 
revision,  but  the  essential  feature,  the  very  small  velocity 
gradients  throughout  the  larger  part  of  the  core  should  remain 
unchanged.  If  it  is  found  desirable  to  change  the  absolute 
travel  times,  this  can  be  done  by  small  velocity  variations 
in  the  inner  and  outer  core  or  in  the  mantle. 

CONCLUSIONS 

The  core  model  presented  in  this  paper  favors  a  division 
of  the  earth's  core  into  four  parts. 

1.  The  inner  core  or  the  DF  zone  between  ■  0-1250  km. 

The  inner  core  has  a  nearly  constant  velocity  throughout  except 
at  the  outer  side  where  a  possible  small  velocity  increase  may 
be  followed  by  a  large  velocity  decrease. 

2.  The  inner  central  region  or  the  BC  zone  between 

"R  1375-2000  km.  This  region  has  a  negative  velocity  gradient 
at  the  inner  end  and  a  125  km  wide  transition  zone  (DC  branch) 
to  bridge  the  gap  to  the  inner  core. 
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3.  The  outer  central  region  or  the  AB  zone  between 
=  2000-2530  km.  In  most  parts  of  this  zone  the  velocity 

gradient  is  small  and  could  be  positive  or  negative. 

4.  The  outer  core  or  SKS  zone  between  ^  *=  2530-3480  km 
The  velocity  gradient  is  assumed  to  be  fairly  constant  and 
somewhat  steeper  than  in  Jeffreys'  model  to  satisfy  the  total 
travel  time  data. 

As  a  final  note  we  should  again  mention  that  the  core 
model  shown  in  Figure  14  is  subject  to  uncertainties  arising 
from  the  problems  of  phase  identifications,  number  of  branches 
and  branch  endpoints.  The  effects  of  the  scatter  in  our  data 
on  the  final  model  were  investigated  (Toksflz,  et  al,  1968). 

It  is  found  that  the  scattering  of  data  does  not  affect  the 
major  features  of  the  above  core  model. 
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Stations  used  for  travel  time  and  dT/dL  measurements 
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TABLE  2 


List  of  earthquakes  used  in  the  core  phase  studies 


Date 

h 

m 

APR 

7,  1964 

13 

18 

MAY 

26, 

1964  * 

10 

59 

SEP 

12, 

1964 

22 

07 

NOV 

1, 

1964 

12 

26 

JAN 

16, 

1965 

11 

32 

JAN 

17, 

1965 

20 

57 

JAN 

24, 

1965* 

00 

11 

MAR 

21, 

1965 

11 

08 

MAY 

16, 

1965 

11 

35 

MAY 

19, 

1965 

06 

03 

MAY 

24, 

1965* 

19 

44 

AUG 

2, 

1965 

13 

19 

SEP 

12, 

1965 

22 

02 

NOV 

21, 

1965 

10 

31 

FEB 

17, 

1966 

11 

48 

Location 


18.9 

0 .  IN 

123. 7E 

12.3 

56. 2S 

2  7.8W 

03.2 

49. IS 

164 .2E 

06.2 

3.  IN 

128. IE 

37.4 

56. 6S 

27. 4W 

41.3 

6.8S 

109. IE 

12.1 

2.4S 

126. 0E 

16 . 2 

1.5S 

126. 5E 

46.0 

5.3N 

125. 7E 

58.9 

6.5S 

105. 4E 

10.9 

56. IS 

27. 6W 

55.9 

56. 2S 

157. 9E 

34.3 

6.3S 

70. 5F. 

54.0 

6.3S 

130. 3E 

00 . 8 

32. 2S 

78. 9E 

Depth 

Mag . 

150  km 

6.30 

120 

6.5 

33 

6.90 

65 

6.3 

101 

6.2 

242 

6.5 

6 

6.6 

33 

6.2 

36 

6.2 

74 

6.3 

120 

6.7 

33 

6.6 

34 

6.2 

33 

6.6 

33 

6.4 

*  No  dT/dA  calculated 
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Figure  Captions 


1.  P  wave  velocity  models  for  the  earth's  core. 

2.  dT/dA  variations  for  the  different  velocity  models. 

2  2 

3.  d  T/dA  plot  for  Adams  and  Randall's  core  model. 

4.  Contour  map  of  travel  time  residuals  from  the  Long  Shot 
explosion. 

5.  Observed  dT/dA  values.  Solid  lines  indicate  the  average 
of  observations  and  used  in  the  velocity  determination. 

6.  Observed  and  calculated  travel  times.  Points  connected 
by  thin  lines  were  interpreted  as  crustal  reverberations. 
Thick  lines  correspond  to  travel  times  calculated  from 
our  core  model. 

7.  Samples  of  P(D'DF)  phases. 

8.  Amplitude  observations  in  the  distance  interval  135°-180#. 

9.  Normalized  amplitude  observations  in  the  distance  interval 
105*-135° . 

10.  Seismograms  showing  P(DC)  phases.  Precursors  are 
indicated  by  the  notation  PRE. 

11.  Sample  seismograms  of  precursor  phases.  Note  the  long 
and  almost  continuous  nature  of  these  waves. 

12.  Samples  of  P(BC)  waves. 

13.  Seismograms  showing  P(AB)  and  P(DF)  phases. 

14.  The  new  core  model.  Velocity  changes  in  the  lowermost 
part  of  the  mantle  are  also  shown. 
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II. 2  Direct  Measurement  of  dT/dA 

(E.S.  Husebye) 

A  very  common  problem  in  seismological  research  is 

the  measurement  of  the  apparent  velocity  or  the  wave 

parameter  p  =  dT/d*~  over  a  seismograph  network.  Several 

methods  here  have  been  developed  (Press,  1954;  Aki,  1963; 

Kelly,  1964  ;  Johnson,  1967;  Knopoff  et  al^.  ,  1967)  which  are 

fully  adequate  for  their  special  purpose.  However,  the 

above  methods  are  limited  in  some  ways,  i.e.  they  are  not 

applicable  to  an  arbitrarily  spaced,  super  large  seismic 

array  (Jansson  and  Husebye,  1967)  where  the  only  known 

parameters  are  relative  arrival  times  and  station  coordinates. 

Recently,  Chinnery  and  Toksttz,  1967,  devised  a  simple 

least  square  method  for  measuring  dT/dA  directly  from  LASA 

data,  assuming  a  plane  wave  front.  Husebye  and  Toksflz, 

1968,  extended  this  method  to  cover  the  general  case  by 

introducing  spherical  coordinates  and  representing  the 

wave  front  as  a  second  order  surface,  i.e.  that  a  Taylor 

expansion  of  T  gives  a  second  order  polynomia  in  4’  (latitude) 

and  (longitude)  whose  coefficients  are  estimated  by 

regression  analysis.  We  then  compute  azimuth,  dT/dA  and 
2  2 

d  T/d^.  .  The  problem  to  be  dealt  with  here  is  to  improve 

the  latter  method  by  transformations  of  the  independent 

variables  $  and  ,i  .  In  short,  our  goal  is  to  express  the 

response  surface  of  T  in  a  simple  form,  and  thus  increase 

the  applicability  of  the  above  method.  It  should  be  noted 

2  2 

that  the  dT/dA  and  d  T/d [  parameters  are  much  more 
sensitive  to  reading  and  identification  errors  than  the 
arrival  time  T. 
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SYMBOLS 


O&l  )  1  x  p  vector  containing  the  unknown  trans¬ 

formation  coefficients 

t<>>  th 

j  elements  of  00;,  for  lb  first  itera- 
^  tion  of  the  jth  transformation 

&  ^  estimated  values  of  06  i 

S  \  a  vector  representing  regression  coefficients 


fi  > 

b  , 

C  , 

&  •> 

£(  )» 
K  i 

*  > 
7 

c-2, ) 
^  , 
T  7 

y , 

1o 

fo 


estimated  values  of  /3 

(k  1-  i).\(k  +  i$  matrix,  c^j  is  an  element  of  C 
epicentral  distance 

unit  vectors  pointing  north  and  east 
expected  vaLues 
expected  value  of  T 
azimuth  angle 


no.  of  observations 


latitude  and  iongitude 
latitude  and  longitude 


variance 


estimate  of  <T 


dependent  variable  or  relative  arrival  time 

1  x  k  vector  containing  independent  variables 

[n  x  (k  +  i)]  matrix  constituted  oi  n  z-vectors 

a  vector  consisting  of  transformed  mements  of  X^ 

a  vector  consisting  of  transformed 
^  and  X 


elements  oi 
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THEORY 


In  the  present  investigation  we  suppose  that  tke  errors 
in  the  T's  are  at  least  approximate Ly  normal  and  independently 
distributed  vltk  constant  variance.  We  skall  concentrate  on 
finding  trans formations  in  the  and  X  t®  reduce  tke  T-function 
to  as  simple  form  as  possible.  Our  aim  is  to  work  wltk  a  poly¬ 
nomial  of  low  degree  in  tke  transformed  variables  ratker  tkan 
a  polynomial  of  klgker  degree  in  tke  original  variable.  Tke 
procedure  outlined  here  is  based  on  a  method  described  by  Box 
and  Tidwell  (1962). 


We  first  start  wltk  a  short  outline  of  Box  and  Tidwell's 
method  which  is  as  follows.  Suppose  observations  Tj/,  i  -  l,...,n 
are  available  at  n  sets  of  conditions  Xl  j  where  Xu. 

is  a  kx  t  vector  giving  tke  levels  of  the  x's  for  the  uth 
observation.  Suppose  further  that: 


E  (Tu)  -  u.  EL(  l  U.  “  *  KvV 


0  i 


(i) 


£ 

where  O'  is  unknown.  Assuming  tkat  tke  reaponse  can  be 
closely  represented  over  tke  region  of  interest  by  some  simple 
function : 


(.li? 


(2) 
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wkere  tke  elemeata  tii*  *  j*fk.  of  tke  vector  ^  are  tke  x's 
treeefenied  la  mm  eultakle  way  ae  r.kat : 


Hi s 


(3) 


iployed,  i.e.  tkry 

-  i*> 


r  l.)  7  .^1 

[•4  VI  >  *  *  '■f>c  J 


wkere  Oi/l  ia  *  p-dimeasioaal  vecter  witk  elemeats  Q^tU 
tke  uakaovra  ceaateata  ef  tke  itk  traas format ioa .  Tke  fuactioaal 
reletioaakip  else  iavelvea  L  uakaewa  coaataata  vkoae  values 
/3-t^  are  tke  elemeats  of  tke  fi  -vector  aad  depead  upoa 

tke  particular  traaa format loaa  of 
depead  oa  tke  ckeice  ef  tke  06'a.  Suppose  <Xf[ 
are  tke  first  fuesses  for  tke  ceastaata  ef  tke  itk  traasfer- 
■atiea.  A  aatural  first  gucsa  would  be  to  use  uatraasforaed 
data,  tkat  is  -  1*0  We  tkea  write  for  tke  vec¬ 
tor  wkose  itk  elemaat  ia  Expaadiag  about 

tkeae  guessed  values  aad  igaeriag  terms  of  aecead  aad  kigker 
orders  ia  (cC ;j  -  ^  ,  we  kave  approximately 

'  J’1  %„:! 


Wo  also  kave: 


W)|HwV 

1 

,  _ 

r  ,  ■) 

_  *3«ui 

*» 

.  ■ 

(5) 


ji;=  ^i0> 


Tke  qusatities 
oquatioa  (3),  wkile  tke  qusatities 


,  are  ebtaiaable  from 

[WVMS’fw1*’ 


must  be  estimated  ia  some  way.  Tki's  caa  be  doae  coaveaieatly 


•  J 
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by  Baking  a  preliminary  fitting  of  the  observations  to 
by  the  least  square  method  or  any  other  suitable  procedure. 
By  differentiating  the  fitted  expression  Tu.’|  ^u.  where 
k  is  the  Least  square  estimate  of  /3  ,  we  may  obtain 

approximate  values  tor  the  quant  ities  (  $  l/V  2  . 

The  n  £  ,  l>'w  approximate  values  of 


( 


o  to 


u) 


Vi 

—  . 

V/ ;  -  w 


(6) 


which  can  now  be  calculated,  provide  a  set  of  "independent 
variables'  fro*  which  we  can  obtain  the  adjustments  to  the 
trans format i on  constants  by  refitting  the  observations  to 
the  whole  expression  on  the  right  side  of  equation  (4).  These 
adjusted  constants  can  now  take  the  place  of  the  "first  quesses 
in  the  above  calculation,  and  the  whole  cycle  is  repeated. 


Scismo  logical  Application 

The  arrival  time  T  of  a  seismic  wave  usually  increases 
monotonica lly  with  epicentral  distance,  and  thus  might  be 
represented  as  a  second  order  polynomia  in  station  coordinates 
^  and  X  .  As  the  earth  is  approximately  spherical,  we 
can  safely  ignore  the  earth  radius  in  the  above  function.  How¬ 
ever,  it  would  be  preferable  to  work  with  an  essential  linear 
function  in  <V  and  X  ,  as  the  number  of  station  records 
available  for  analysis  usuaLly  is  limited,  say  5-15.  Note 
that  the  above  approach  makes  us  quite  free  to  consider  any 
number  of  seismologies  1  stations  as  an  array  which  we  might 
use  for  measuring  dT/dA  .  As  array  geometry,  azimuth,  epi¬ 
central  distance  and  the  wave  under  consideration  will  vary, 
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tke  function  which  best  fits  tke  observed  response  surface 
T  will  vary  too.  Tke  Box  aad  Tidwell  method  outlined  above 
should  therefore  by  appropriate  for  direct  Measurement  of 
dT/d  h  ,  as  the  response  surface  kere  is  not  tied  to  a  single 
class  of  functions.  For  example 


L.X.  **  ‘  ° 


(7) 


includes  many  forms  commonly  found  useful,  like  square  root 
and  reciprocal  transformations.  Tke  above  procedure  can 
easily  be  adapted  to  transformations  of  tke  dependent  variable 
too.  For  example,  assume  we  make  guesses  different  from  1.0 
for  tke  values  of  06i,j  in  tke  first  iteration.  We  can  always 
regard  these  new  variables,  Including  any  transformation  of 
tke  dependent  variable,  as  being  tke  basic  variables.  Wken 
a  second  iteration  is  performed,  only  tke  Xi  *  will  be  taken 
as  new  basic  variables. 

Leaving  the  general  notation  above,  and  returning  to 
our  problem,  namely  to  find  the  response  surface  of  the  rela¬ 
tive  arrival  time  over  an  arbitrary  array,  we  then  have: 


\  =  fio  Jo  y3|  To  ‘  Z  k. 


(8) 


where  Z0  is  a  dunmy  variable  always  equal  l  •  0^  Z  \ s  ^  ** 

etc.,  k  is  number  of  parameters  in  the  T  expansion.  Our  first 
guess  is  that  z  *  /•  o  ,  and  rewriting  of  eq .  (4) 
gives 


fu..  i  +.  (ot,  -  to)  Sr/fo,  e>> /«w 


(9) 
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Fro«  eq  .  (5)  we  have: 

7^-T  +  I o)  ^  Z&y  •Xb'X  (10) 

We  start  by  applying  the  Method  of  least  square  to  eq .  (8). 
Differentiating  this  equation  with  ^respect  to  41  and  A  , 
we  May  calculate  *9 T / c )(^  and  9t/9\  The  sign 

roof  aeans  statistical  averaging.  Construction  of  new, 
independent  variables  gives: 

zk+-t  ^VsV^'k'' 'V  UO 

When  refitting,  we  use  the  equation: 


lu.  -  /io^o  +  •  •  ■+  (061  -  \  +  ^12) 

Thus  we  get  esti*ates  and  '  l*o)  of 

and  (ftLx-l-o^)  .  The  r  and  nay  be  treated  as  ^ 


and 


\ 


in  further  iterations. 


ADEQUACY  OF  TRANSFORMATION 


It  is  desirable  to  have  sone  Measure  of  precision  of  the 
(X  t  quantities,  both  for  deciding  nuMber  of  iterations  and 
whether  it  is  worthwhile  to  transfom  the  original  variables. 
The  variance  of  CU  jV<*h  (<Xi)  might  be  calculated  in  the  ordi¬ 
nary  way  for  estinated  least  square  coefficients.  Laagining 
Matrix  notation  applied  to  eq .  (12),  then 

Vdh(.CU)  -  (13) 


where  0^1,^ 


the  matr  *  (y'y)  »  is  the  variance  of  T,  and  its  esti- 


is  the 


(lou) 


diagonal  element  in 


lated  value  is  the  residual  Mean  square  value  s  of 


A 

T  • 
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However ,  when  the  number  of  observations  is  snail, 
allowances  are  necessary  for  the  -astimates  Included  in 
the  variables.  When  working  with  transformed 

linear  functions^  this  difficulty  may  be  avoided  by  using 
Fieller's  theorem  (1940).  However,  the  straight  forward 
method  here  is  to  test  the  transformations  by  calculating 
standard  parameters  like  sum  of  squares  around  regression 
and  due  to  regression,  lack  of  fit,  etc. 

CALCULATION  OF  AZIMUTH^  dT/d  A,  AND  d2T/d&  2 


When  the  above  steps  are  finished,  the  travel  time 
variation  over  the  array  is  expressible  as: 

T  *  l0Zo  4  •  •  ■  ■  •  H  ^  ^ 

where  Ze*  Lo)  Zi  s 
variance  of  is,  setting  c  -  [y 

VarCu)  =  C%  (V  C  i:)  U5) 

wher^  "Z c.  =  ^  ^  ^  •  •  ••  2lc^). 

Differentiating  the  scalar  T^,  remembering  that  we  are  using 
spherical  coordinates,  we  get  if  is  a  second  order  poly- 
nomla : 


^Vie^  t 


(16) 
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First  order  expansion  of  (16)  gives  (  LL  -tern  only): 

Ul*»  dUL  *|l>i4cU>i+  3.(byvJ^>)^fc  +  (17) 

Since  we  have  assuned  that  the  error  ATi,  accompanying  the 
ith  observation,  is  an  independent  Gaussian  variable,  having 
zero  nean  and  variance  0'^' ,  dLw  and  da  i  will  also  be 
Gaussian  variables  with  zero  neans.  Performing  statistical 
averaging  over  eq .  (17),  we  have  for  the  variance  of  Hi.  : 


Va>-  (Ui)  di, .  $ 


-  l  sOj 


(18) 


Since  we  are  averaging  over  the 


dt>i  t 


erms ,  we  nay 


write : 


V^Ou)  ~|Vair  (t,)  +  h) 

+  H Sp  Cov  ( U>)  ♦  i X^Ccv  (  U  U) 

T  Cov(ti^U)  j 


(19) 


where 


Cov  (,  V'i) 


means  covariance  t\i^3  and  corresponds 


to  <j  Multiplied  with  the  C\-^  elenent  of  the  C-natrix 
given  above.  In  similar  ways  we  nay  calculate  the  terns 

„d  C.v  (lUV;.)  . 


The  azimuth  of  the  wavefront  is: 

W  yt 


_Vi_ 

UU 


(70) 
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First  order  expansion  of  both  sides  of  eq .  (20)  gives: 


tanVi.  ,  1  h.  +  ikijy.1'  ~  jfc ^  — 

*  <**  Uc 

Using  eq.  (20)  we  get: 

aVu  4JU 

Lv.  aw 


(21) 


x  Cos'jfi  rSih^Co^C, 


(22) 


(23) 


Thus  the  error  in  ^ o  is  Gaussian  with  zero  scan.  The  vari 
ance  of  ^ *  is  : 

Vd H  (  Vi)  =  _ ICoV  ..Vat- Pi) 

‘ 1  1  1  rw5^  ucv.  ^jar 

The  wave  parameter  -  at/al  is  : 

:p;N  LL1'  *  V;* 

Differentiating  eq .  (24),  we  get: 

UidULi  +  (25) 


(24) 


Als 


(26) 


o  al  ^  i  is  Gaussian,  and  the  variance  of  ^  is: 

Va>l?0  =  ^ 1  [ill1  VarLUi) U  Cov(lUV^V^Va.^j 

The  dfV&k  or  term  is  very  important  in  interpreta¬ 

tion  of  amplitude  data,  and  it  might  be  worthwhile  to  calculate 
when  our  observations  are  very  precise.  Normally,  this  quan¬ 
tity  Is  obtained  from  dT/d  L  observations,  but  for  the  sake 
of  completeness ,  the  calculation  procedure  will  be  outlined. 

By  taking  the  directional  derivative  of  the  p-vector,  we  get 
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2  ? 

for  d  T/d  tk  -  vector: 


Jt-tyf  -  [L- ?]  ^ e**  • 

-  e«r-%si§A  ■[llic*  +  V,iSxl 

"  co5^9\  J  L  t  J 

Introducing,  new  vectors  L  and  M:  - 

Proceeding  in  the  sane  ir.anner  as  we  did  foi  calculating  the 
^—vector,  we  have: 

Az.n(^T/cJtf)-  A^Ulrv  Ov/u) 


(27) 


(28) 


(29) 


Vd  +  (Aili)  r  Si^  /Uv  tCo  £  f\ZiL 


j/dV-CHc)  _  ZCovCUHl)  +Va<r  (IQ 

L  ft?  Me  L  Ui7, 


ana  ,  r  ,1  iy 

V^^’W'  [L1  VWu)+  Z  UH^Cov  (I4^VM^ 


(30) 


A  large  problem  in  present  day  seismogram  analysis  is 
the  existence  of  local  travel  tine  anomalies.  This  effect 
should  be  properl>  accounted  for  to  avoid  biased  errors  in 
the  dT/dA  estimates. 

If  the  epicenter  parameters  are  not  given,  these  could 
be  easily  calculated  from  the  estimated  azimuth  and  dT/dA 
values  if  a  dT/d  A  table  is  available.  This  has  been  demon¬ 
strated  by  Kelly  (1964)  and  the  procedure  here  is  as  follows. 
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We  construct  a  spherical  triangle  whose  corners  and  corner 
angles  ara  $  T  A(  /lU ;  jfO  i 

and  hi? OLR^O.}-  •  The  angular  lengths  of 

the  sides  are  From  these  quantities 

and  the  following  equations,  we  can  easily  calculat  ^ 
ai.d  A  . 

Co^(Ai'A)  -[co5<f<: COSMOS Ai t0^l/ 

Due  to  existence  of  biased  azimuth  errors,  it  night  sometimes 
he  preferable  to  combine  the  above  method  with  another  one 
described  by  Husebye  (1966).  The  latter  method  needs  some 
modification  for  analytical  usage. 

Demonstration  of  the  Method 


The  above  method  has  been  applied  to  the  two  arrays 
shown  in  Fig.  1.  The  P  travel  time  data  have  been  taken  from 
Jeffreys -Bullen  tables  and  in  some  cases  random  errors  have 
been  added.  Two  types  of  response  surfaces  are  investigated. 


T-k  + 


(32) 


7=  L.A.f  +Mai*Uf^+WfM+  <33> 

Input  data  like  station  coordinates,  travel  times  and  also 
true  dT/d&  and  azimuth  values,  are  given  In  Table 
The  results  are  displayed  in  Tables  2  and  3,  and  need  a 
few  comments.  For  the  second  order  response  surface  given  hy 
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eq .  (33),  no  improvement  is  gained  by  the  iteration  process. 
The  new  Ob  -coef ficents  are  very  close  to  1.0  having  stan¬ 
dard  errors  around  0.6  .  This  result  reflects  the  very  nature 
of  our  problem,  i.e.  to  measure  the  slope  of  a  surface  which 
deviates  slightly  from  a  shell.  For  the  first  order  response 
surface,  the  new  Ob  -coef ficents  give  sometimes  a  small  im- 
proveraent  in  the  dT/d&  values.  In  these  cases  it  was  suf¬ 
ficient  with  only  one  iteration.  Implemented  in  least  square 
analysis  is  that  the  most  accurate  estimate  of  the  sought  para 
meters  corresponds  to  points  at  the  center  of  the  array,  so 
only  these  values  should  be  retained  for  further  investiga¬ 
tions.  Important,  the  "center1  values  in  Tables  2  and  3  ex¬ 
hibit  little  dependence  both  on  the  type  of  response  surface 
and  inclusion  of  random  errors.  For  example,  analysis  of 
core  phases  could  sometimes  be  performed,  using  only  two  free 
variables  in  the  response  surface. 

Discussion 


The  method  outlined  above,  favors  a  somewhat  different 
approach  to  the  analysis  of  travel  time  data.  Instead  of 
using  many  events  recorded  at  widely  scattered  stations,  it 
is  preferable  to  use  fewer  events  recorded  in  a  relatively 
dense  station  network.  If  digitized  data  are  available,  the 
proper  procedure  might  be  as  follows.  The  first  step  is  to 
band  pass  and  remode  filter  the  records  if  this  is  deemed 
necessary.  To  avoid  the  non-linear  properties  of  the  general 
remode  filter  (Flinn,  1965),  an  acceptable  version  of  this 
would  be  to  rotate  the  seismograph  system  such  that  one  of 
the  component  points  in  the  direction  of  the  particle  motion 
of  the  phase  under  investigation.  By  weighting,  the  ellipti- 


-52- 


city  of  the  particle  notions  is  restricted  to  certain  Units  , 
and  parts  of  the  reverberation  effects  would  be  suppressed. 
Th'i  first  order  tine  lags  are  calculated  fro*  proper  know¬ 
ledge  of  focal  parameters,  while  second  order  tine  lags  are 
obtained  by  using  a  cross-correlation-iteration  scheme  (Fair¬ 
born  and  Gangi,  1967,  see  aLso  Jansson  and  Husebye,  1968). 
Additional  refiennent  of  the  arrival  tine  values  are  possible 
when  information  on  local  velocity  anomalies  or  station  cor¬ 
rections  are  available.  Essentially,  the  above  procedure 
corresponds  to  the  much  used  method  of  velocity  filtering  of 
data  from  small  arrays. 
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Figure  Captions 


Fig. 1.  Array  configurations.  The  triangles  indicate 
the  centers  of  the  arrays. 


Table  Captions 

Table  1.  Input  data  used  in  the  application  and  testing 
of  this  method. 

Table  2a.  Calculated  azimuth  angles,  dT/dA  values  and  time 
residuals,  assumint  an  eq.  (32)  response  surface 
model. 

Table  2b.  Calculated  azimuth  angles,  dT/dA  values  and  time 
residuals,  assuming  an  eq.  (32)  response  surface 
mode  1 . 

Table  3.  Calculated  azimuth  angles,  dT/dA  values  and  time 
residuals,  assuming  an  eq.  (33)  response  surface 
mode  1 . 


Sta.  Coordinates  Distance  Azimuth  Travel  dT/dA  Model  for 

no.  North  East  (deg)  (deg)  time  (sec/deg)  random  errors 

(deg)  (deg)  (sec)  in  travel  time  (sec) 
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III.  THEORETICAL  SEISMOGRAMS  AND  EARTH  STRUCTURE 

III.l  Upper  Mantle  Seismic  Models 

(D.V.  Helmberger) 

In  this  study  we  will  generate  synthetic  seismograms 
for  upper  mantle  models  of  spherical,  homogeneous  layers. 
The  seismic  signals  were  generated  by  an  underground  shot 
at  the  NTS.  We  will  synthesize  the  motion  at  the  variour 
LRSM  stations. 

The  procedure  followed  will  be  to  develop  the  theory 
for  a  plane  layered  earth  assuming  a  point  source.  Next, 
the  theory  of  spherical  layers  will  be  treated  and  shown 
to  reduce  to  a  simple  modification  of  plane  layered  theory 
when  the  curvature  is  small.  Finally  this  theory  is 
applied  to  models  proposed  from  observations  of  the  Bilby 
event.  A  short  comparison  between  these  observations  and 
synthetic  waveforms  is  given. 
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Point  source  in  s  plane- layered  earth 


Suppose  we  consider  the  problem  setup  specified  in  the 


following  die  gran.  ^ 

p  •  *  f  t  5  « ►*«'  c  ^ 


,^ct  l^tr 


ta,  c,  ?,  P, 


_  _  ______  rj,*_  _cn  s„  n,. 

Tk  C  ~  P 

M  f  /  t  *  *  »  f  l 

Let  W  be  tbs  vertical  displacement  with  r  the  horizontal 
coordinate.  The  Laplace  transfora  with  respect  to  tiae  is 
denoted  by  an  overbar  with  a  the  transfora  variable.  Sup¬ 
pose  the  point  source  emitts  a  step  function  of  unit  strength, 
Then  the  integral  expression  for  the  transformed  response 
from  the  bottom  of  the  nth  layer  is 

—  (  **  —  "s  3<r>  „  (fi,  I ) 

wli-.s)  s  cju  kc(*t*)#tr)#ir)  *(r>  e  ^  J 


where 


<1  <f>  =  ^  r*-  %■ 


ri,  * 


V(r‘ 


r:  ■  r*> 


I 


&<f)  -  Reflection  Coefficient  at  n,n+l 

j  j  ■*  j  j 

The  above  formula  yields  the  transformed  displacement  for 
the  generalized  ray  that  has  traversed  the  strata  exclusively 
in  the  P-mode  and  contains  no  internal  reflections.  The 
generalized  ray  that  is  reflected  at  the  surface  above  the 
source  will  be  included  in  the  source  function  and  will  be 
discussed  later.  T(p)  and  R(p)  are  the  transmission  and 
reflection  coefficients  and  are  given  by  Helaberger  (1968). 
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The  formalism  applied  In  this  section  is  obtained  directly 

from  that  paper.  Equation  (Al)  should  be  compared  with  (1), 
The  extra  s  and  D(p)  are  introduced  by  the  fact  that  the 
receiver  is  situated  at  the  surface  and  because  we  are  in¬ 
terested  in  displacement,  see  Knopoff,  Gilbert,  and  Pliant 
(1960)  for  details. 

Applying  the  inverse  Laplace  transform  we  obtain 

t 

w(r,t)  =  2.  A.  r(lsL\e<£Wj^JiL*zlLLL 

*  jt  J  '  tr I, _  _  , -j /  .  l t 

where  (t  -  *  >  (  t  -  T  +  »  f  ►  ) 


Tif)  =  +  a  ^ 

f/i)T  =  (  r  -  a  f  r  V/T  ) 


The  j's  are  to  be  summed  from  1  to  n.  Suppose  we  now  apply 
the  high  frequency  approximation  which  is  valid  for  sources 


of  short  duration 


Then  (A2)  reduces  to 


%  -  ?  +  7  p  t  3/ 


W  (\r,  t  )  = 


=  *.  i-f  *  +  ft  >  1 

ir  j  t  L  \Tt^  J 


where 


*<t)  =  oU  )~==  ) 

To  compute  the  vertical  displacement  we  must  sum  over  all 
the  generalized  rays  which  arrive  in  soHe  specified  time  or 


where 


(  r,  t  )  =  / 

=  7  A  (rf  * 
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..  .  -re  numbered  1  to  n  see  Helmberger 

where  the  generalized  ray.  are  nuaberen 

(1968)  for  detail*.  Now  one  can  perfor.  the  au—tlon  fir.t 
yielding 


Wr  (f,  t )  =  £  £  [  re  *  4  f” (  *  ' ] 

<  _  1 _ 


(A*) 


The  high-frequency  approximate  solution  ha.  been  compare 
With  the  exact  solution  and  found  very  .ati.factory  for 
body  wave  analysis  when  the  source  duration  is  snort  com¬ 
pared  with  the  travel  time  between  the  source  and  receiver. 
This  approximation  allows  o*  to  replace  n  convolutions  by 
one.  Before  generating  synthetic  seismogram,  we  must  con¬ 
sider  the  effect  of  curvature  which  we  now  take  up. 


B.  Point  source  in  a  spherically  layered  earth 


We  first  consider  the  problem  solved  by  Gilbert  (1960) 
on  scattering  of  impulsive  SH-wave.  fro.  a  rigid  sphere. 
The  problem  setup  is  tfiven  below. 


The  displace  ent  *4 
coordinates  (r,  ©* $ 


has  only  a  cotaponer.c  in  spherical 
).  Suppose  we  assume  a  delta  function 


ti.  e  depend.; ul  source,  then 


j * u  -  -  S(r-r0) 


(8/) 


v  u  - 


y3x  J  t  x 


y  *  0 
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wliere  A  is  the  velocity.  If  we  apply  the  Laplace  trans¬ 
form  on  the  rime  variable  and  the  Legendre  transform  on 
the  angle  variable  o  we  obtain 


J-  ±  f  r'ii'j  -  ( 

r*  d?  r  L  <9  r\  \  A 


_L  J(J-ti))  U  =  -_J_  J(y-y  ;  {BS ) 
Y  1  J  ^  * 


where 


V 


(v)  =  [  u(c~*o) 


and 


ot 


a  i  )  * 


(b*; 


where  the  L  contour  runs  above  the  real  /?  axis  frou 
to  <*>  .  A  particular  solution  of  (B2)  is 


C/o  =  *  y  )  ^  ( k  r0  )  r  <  j  *t  -  S'A 

where  and  kr  are  the  spherical  bessel  functions.  Since 
the  sphere  is  rigid  l/=o  at  k-  a  and  the  total  solution 
becomes 


U=  -  l±  )  krUrr) 

'iff*'  K  - 

*,  (*‘> 

We  will  be  primarily  interested  in  the  second  term  which  re¬ 
presents  the  reflected  motion.  The  solution  in  the  above 
form  is  due  to  Gilbert  (1960).  He  applies  the  saddle  point 
approximation  to  obtain  the  first -motion  solution.  However, 
we  can  obtain  the  high  frequency  approximation  by  the  appli¬ 
cation  ofCagniard's  method. 

Let 

«.  tr  i  , 

P  ,  C  ot*>  e  )  —  «.  tu  PC  c*~t* )  —  3  «.  Ur^’  *  *  Q  (  a ) 

A'1 


and  ci  =•  tr,  +  m,  t  where 


u 


l 


(  CO  *  t?  ) 


u(J-i)** 
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and 


az  =  -  [  ^  o' (c~*;  v(j-  jt)  ^  f-; 


We  will  consider  k.*  which  yields  the  reflected  notion. 

,  -<  {IV-  % ; 

_ _  .  J  ->  ^ 


V,  =  - 


and  let 
J  =  Lm  +  Q  = 
Equation  (B4)  becomes 


(  a  Jr^  S/n  $J  1 


-  ^ 


0 


I 

C 


/<r  2  $  J 


6<z  C  6  ) 


_  _1  f 

f  * 


-y*4  e 


c  c  (  a  fr si  «e  J 


'h 


.  i 


where 
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The  asyicptotic  forms  of  the  bessel  functions  were  used  in 
evaluating  V  .  We  also  used  the  following  identity 
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Substituting  we  obtain 
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Now  let  ^  *  s  p  and  change  variable,  (B5)  becowes 
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and  we  are  now  assuming  a  step  function  source.  Let  ~p>  1  =  /■//" 
and  (B6)  becomes 
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where 
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Now  applying  Che  Laplace  inversion  we  note  that 
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because  x  «  rs~n©  .  The  form  of  (B7)  is  identical 
to  (A3)  although  (  /*t  )  is  different. 

This  solution  can  be  generalized  to  n -layers  when  we  neglect 
terms  involving  (as)1 
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normalized  radius  of:  layer  -j  and  f  is 

replaced  by  (f  >  )  in  all  other  functions  of  the  plane- 
layer  case  functions  where  *  is  the  normalized  local  radius. 
That  is  at  interface  K  the  transmission  coefficient  is  the 
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same  as  that  in  the  plane  case  except  f  is  replaced  by  tftr. 
The  solution  is  only  rigorous  for  reflected  responses  al¬ 
though  it  is  probably  good  for  refractions  as  well  when 
the  curvature  is  saiall. 


C.  Synthetic  seiSB:ograns  based  on  the  NTS  Bilby  event 

Synt  etic  seist.ograais  are  produced  by  a  double  convolu¬ 
tion 

SS  (  t  )  -  JL  I"  Wit)  *  T(t)J 

where  _ 

T(t  )  =  [  S it )  *  T(t)  J 

and  X(tj  is  the  instrumentation  response  and  S It)  is  the 
idealized  source  function.  We  chose  to  include  the  inter¬ 
action  of  the  earth's  surface  with  the  direct  wave  and  label 
the  total  radietion  S (*).  At  ranges  greater  than  30  degrees 
the  observations  are  quite  similar  in  character  which  reflects 
the  nature  of  Tit)-  This  is  expected  because  the  model  struc¬ 
ture  is  apparently  simple  at  depths  greater  than  900  km  and  Wit) 
is  approximately  a  step  function  so  that  SS(t)  s:  Tit) 

The  point  will  become  clearer  when  we  cov  pute  a  simple  model 
(S.E.  Profile).  We  have  done  some  calculations  on  Tit)  directly 
from  XU)  and  Sit/  allowing  *  layer  above  the  source. 

Using  the  X(tJ  discussed  by  Carpenter  (1967)  and  source  function 
for  the  Bilby  event  by  Toksbz  (1967),  we  calculated  T*  dis¬ 
played  in  Fig.  1  where  we  assumed  a  thin  layer  of  low  velocity 
material  above  the  shot.  The  theoretical  response  is  the  mo¬ 
tion  appropriate  for  energy  emerging  at  about  30  degrees.  The 
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decaying  oscillation  effect  is  caused  by  energy  leaking  out 
of  the  layer.  The  lower  plot  is  Tit)  which  contains  a  con¬ 
volution  of  the  upper  plot.  Tit)  does  not  change  with  angle 
very  rapidly  so  that  Tit)  is  appropriate  for  a  10  degree 
range  as  well.  Fig.  2  displays  a  profile  of  observations 
running  S.E.  from  the  NTS  taken  from  the  Bilby  event.  Their 
locations  are  given  in  Fig.  3.  Note  that  the  first  3  peaks 
of  the  Oft  observation  look  similar  to  Fig.  1.  However, 
we  do  not  know  Sit)  very  accurately  so  that  we  will  use  the 
first  three  secs  of  the  observation  at  OR  as  the  transfer 
function  Tit)  in  this  preliminary  study.  We  will  discuss 
the  N.E.  profile  first  and  return  to  the  S.E.  later. 

The  N.E.  Profile  has  been  studied  by  Julian  and  Anderson 
(1963).  They  find  a  model  which  is  quite  similar  to  John¬ 
son's  (1967).  A  simplified  model  which  contains  the  most 
prominent  features  is  given  in  Fig.  4.  Synthetic  seismo¬ 
grams  for  this  model  are  displayed  in  Fig.  5  and  should 
be  compared  with  the  observations  in  Fig.  6.  It  should 
be  noted  that  not  all  of  the  observations  are  from  Bilby 
and  difference  in  transfer  functions  are  expected.  Never¬ 
theless  the  prominent  second  arrival  occuring  at  about  12 
degrees  is  quite  impressive.  This  arrival  is  the  critical 
reflection  off  the  400  km  transition.  At  greater  ranges 
this  arrival  becomes  the  first  event  followed  closely  and 
interfering  with  the  response  from  the  b50  km  transition. 

This  composite  waveform  is  also  indicated  in  the  observa¬ 
tions.  At  the  smaller  ranges  the  model  produces  a  rather 
large  arrival  which  is  the  sum  of  reflections  off  the  bot¬ 
tom  of  the  low  velocity  zone.  The  observations  at  this 
range  (not  shown  in  Fig.  6)  do  not  show  this  event.  This 
means  the  model  needs  a  smoother  transition  at  this  depth. 
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We  used  the  model  proposed  by  Dowling  and  Nuttli  (1964) 
for  the  S.E.  profile.  Their  model  la  given  in  Figure  7. 

We  used  the  Birch  model  for  the  shear  velocities  and  den¬ 
sities.  The  velocity  depth  function  was  approximated  by  a 
stack  of  layers  20  km  thick.  Only  the  P-mode  generalized 
rays  were  allowed.  The  rays  containing  only  one  reflection 
dominate  the  notion  although  internal  reflections  were  con¬ 
sidered  near  critical  angles.  That  is  we  include  internal 
reflections  from  five  layers  above  and  below  that  depth 
where  the  geometrical  ray  bottoms.  A  profile  of  synthe¬ 
tic  waveforms  are  compared  with  corresponding  observations 
in  Fig.  8.  The  £S(t)  display  the  fact  that  linear  gradi¬ 
ents  return  the  Incident  waveform  with  very  little  distor¬ 
tion  of  shape,  as  is  well  known.  The  rate  of  decay  as  a 
function  of  range  is  also  quite  smooth  although  there  is 
a  change  of  slope  in  the  velocity  depth  function  at  about 
500  km.  The  transfer  function  is  quite  apparent  in  the 
observations.  However,  the  amplitudes  disagree,  alse  the 
2274  km  observation  has  a  definite  precursor.  This  indi¬ 
cates  that  the  650  ku  transition  is  also  present  in  this 
profile.  The  critical  reflection  occurs  at  about  this  range 
which  would  explain  the  large  amplitude.  None  of  the  S.E. 
recordings  show  prominent  second  arrivals  so  that  the  400 
km  transition  must  not  be  as  sharp  as  in  the  N.E.  profile. 
Thus  it  appears  that  the  upper  mantle  is  quite  inhomogeneous. 

We  are  now  in  the  process  of  gathering  observations  from 
many  shots  in  an  attempt  to  give  a  detailed  treatment  of 
these  transitions. 
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figure  1  Upper  plot  is  the  response  of  a  point  eource 
5  interacting  with  a  surface  layer  assuning  a 

•cep  function  source.  The  lower  la  the  result¬ 
ing  Tit) 


Figure  2  Profile  of  observations  fro*  Bilby 

figure  3  location  of  Bilby  stations  (after  Dowling  and 
Nutt1 i) 


figure  4  Medoi  of  the  upper  nantle 

figure  5  Synthetic  selsaegrans  based  on  the  nedel  in 
fig.  4 

figure  6  Record  section  for  NTS-NE  profile  (after  Julian 
and  Anderson) 


figure  7  Velocity  eiodel  for  data  SB  frow  Bilby  (after 
Dowling  and  Nuttli) 

figure  8  Comparison  of  synthetic  and  observations 
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A  Travel  Time  and  Amplitude  Interpretation  of 
a  Marine  Refraction  Profile 

(D.V.  Helmberger  and  G.B.  Morris*) 

III. 2-1:  Primary  Waves 

In  recent  years  with  the  advent  of  modern  computers 
it  has  become  possible  to  gain  meaningful  information  from 
seismic  amplitudes.  Some  of  the  difficulties  in  this  work 
have  been  obtaining  adequate  descriptions  of  source 
functions  and  maintaining  a  favorable  signal-to-noise  ratio 
at  the  receiver.  It  appears  oceanic  investigations  have 
an  advantage  in  the  first  of  these  whereas  continental  work 
has  the  advantage  in  the  latter.  However,  last  year  at 
the  proposed  Mohole  site  off  Hawaii,  both  features  were 
well  satisfied.  Observations  were  recorded  by  FLIP,  the 
Scripps  Institution  of  Oceanography's  Floating  Instrument 
Platform,  which  provided  an  extremely  stable  platform  with 
a  low  noise  level.  Seismic  signals  were  produced  by  85  kg 
charges  fired  at  roughly  constant  depths  along  a  120  km 
line.  The  observed  set  of  recordings  provided  an  excellent 
opportunity  to  study  the  amplitude  of  seismic  head  waves. 

The  general  method  of  study  consisted  of  first 
obtaining  a  simple  layered  model  which  fit  the  travel  times 
of  the  observed  records.  Theoretical  synthetic  records 
were  computed  for  this  model,  after  which  waveforms  and 
amplitude  decay  rates  of  the  synthetic  and  observed  records 
were  compared.  Where  there  was  marked  disagreement, 
adjustments  were  made  to  the  original  layered  model  by 
trial  and  error.  This  adjusting  was  continued  until  the 
synthetic  records  displayed  the  same  pertinent  features  as 
the  observed  records.  Such  a  technique  is  subjective  and 
does  not  necessarily  yield  a  unique  solution.  Nevertheless, 
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by  fitting  both  amplitude  and  travel  time  the  solution  is 
more  restricted  than  a  simple  travel  time  interpretation. 

OBSERVED  RECORDS 

The  refraction  records  used  in  this  study  were  taken 
during  the  summer  of  1966  north  of  the  Hawaiian  Islands. 

The  particular  set  of  records,  station  SH31,  located  350 
km  north-northwest  of  the  island  of  Hawaii,  is  only  one 
of  the  many  profiles  taken  by  a  joint  operation  of  Scripps 
Institution  of  Oceanography,  University  of  Hawaii,  Oregon 
State  University,  and  the  University  of  Wisconsin.  Aboard 
FLIP  the  outputs  of  three  hydrophones  were  recorded  with 
the  usual  instrumental  set-up  used  by  Scripps  and  described 
by  Raitt  (1956)  and  Shor  (1963).  A  sample  recording  is 
given  in  Figure  1.  The  first  arrival  Gx  is  identified  as 
a  primary  wave.  The  second  arrival  G ^  arrives  with  a 
velocity  of  3.5  km/sec  and  is  tentatively  identified  as  a 
converted  shear  wave;  that  is  a  wave  that  travels  through 
the  ocean  as  a  primary  wave  and  through  the  solid  crust  as 
a  shear  wave.  Figure  2  shows  all  the  data  points  together 
with  a  layered  model  which  is  based  on  the  time  information. 
All  of  the  shots  regardless  of  charge  size  were  used  for 
this  initial  interpretation,  however,  only  the  85  kg  shots 
were  used  in  the  later  comparisons.  We  should  mention 
that  the  layered  model  was  not  based  solely  on  the  information 
from  this  one  profile.  This  profile  was  one  of  a  set 
which  enabled  us  to  compute  split  and  reversed  solutions. 

The  layered  model  shown  in  Figure  2  is  a  solution  which  is 
compatible  with  this  additional  time  information. 

AMPLITUDE  MEASUREMENTS  OF  OBSERVED  RECORDS 

To  compare  the  amplitude~distance  decay  rates  of  the 
synthetic  and  observed  records,  we  had  to  arrive  at  some 
measure  of  the  amplitude  of  the  received  events.  These 
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amplitudes  were  measured  only  on  the  records  of  the  85  kg 
charges.  The  records  for  smaller  charges  were  not  used  in 
order  to  avoid  empirical  corrections  for  charge  size  and 
depth  of  detonation.  The  pulse  amplitude  was  chosen  as  the 
average  of  the  heights  from  the  first  valley  to  the  second 
peak  and  from  the  second  valley  to  the  third  peak  for  the 
three  hydrophones.  The  first  pressure  peak  was  ignored 
as  it  was  usually  too  small  to  measure  accurately.  These 
averaged  values  after  being  adjusted  for  amplifier  gain 
setting  are  plotted  for  two  seismic  lines  in  Figures  3  and 
4.  Figure  3  is  for  profile  SH31,  whereas  Figure  4  is 
for  SHll,  a  profile  some  450  km  northwest  of  SH31.  Figure 
4  is  included  to  support  the  premise  that  other  lines  in 
the  area  exhibit  the  same  general  amplitude-distance  decay 
relation  as  SH31. 

In  making  the  amplitude  measurements  we  observed  that 
the  principle  cause  of  scatter  seems  to  be  interference 
between  the  various  waves.  Most  of  the  events  whose 
amplitude  appears  anomalously  strong  or  weak  on  the  plots 
also  have  a  "disturbed"  waveform  suggesting  that  they  are 
the  result  of  more  than  one  principle  reflected  or 
refracted  event.  This  observation  is  also  supported  by 
noticing  on  the  amplitude  plots  that  the  strong  or  weak 
second  arrivals  do  not  coincide  with  the  strong  or  weak 
first  arrivals,  which  would  be  the  case  if  an  entire  record 
was  anomalously  strong  or  weak. 

Regardless  of  the  scatter  of  points,  the  amplitude- 
distance  curves  display  several  prominent  features.  Out 
to  distances  of  about  30  km,  which  is  the  crust-mantle 
crossover,  the  amplitude  decays  slowly  with  distance. 

Beyond  this  cross-over  point  and  out  to  about  50  km  the 
amplitude  shows  little  or  no  significant  decay  with  distance. 
An  additional  feature  not  shown  on  the  amplitude-distance 


plots  is  exemplified  in  Figure  1.  This  record,  taken  at 
18.8  km  from  the  source,  shows  refracted  and  reflected  events 
coming  from  the  upper  layers.  The  pertinent  feature  is 
how  the  energy  is  distributed  with  time.  The  energy  in  the 
first  arrival  G^  occurs  in  the  first  three  cycles  after 
which  there  is  a  low  signal  level  followed  by  a  prominent 
second  arrival  G2.  It  is  this  distribution  of  energy 
associated  with  the  upper  layers,  together  with  the 
amplitude-distance  relation  which  we  will  attempt  to 
explain  with  the  theoretically  computed  synthetic  records. 

SYNTHETIC  WAVEFORMS  AND  COMPARISONS 

The  procedure  for  computing  the  synthetic  records  is 
similar  to  that  developed  by  Helmberger  (1968).  In  that 
paper  the  synthetic  response  was  obtained  by  convolving  the 
unit  step  function  response  of  the  earth  model  with  the 
derivative  of  the  transfer  function  with  respect  to  time. 

The  model  response  was  calculated  for  a  point  source  and 
receiver  in  the  oceanic  liquid  layer  overlying  a  crustal 
model  of  plane,  homogeneous,  isotropic  solid  layers  and 
a  solid  half  space.  The  method  of  generalized  reflection 
and  transmission  coefficients  was  used  to  calculate  the 
transient  response  in  integral  form.  Such  integrals  can 
be  solved  by  the  Cagniard  (1939)  method.  A  version 
modified  by  deHoop  (1960)  and,  later,  extended  to  point 
sources  by  Gilbert  (1963)  was  used.  The  transfer  function 
or  a  modified  source  function  which  includes  the  effects  of 
bubble  pulses  and  various  reflections  from  the  surface  was 
defined  as:  . 

T(t)  -  W(t)  -  W(t  -  tx)  -  W(t  -  t1  -  t2) 
where  W(t)  is  the  convolution  of  the  source  function  and 
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the  instrumentation  response,  t^  and  t2  are  time  lags 
determined  by  the  depth  of  source  and  receiver  respectively. 
The  method  of  analysis  used  here  is  the  same  except  for  the 
addition  of  an  absorption  factor. 

In  order  to  avoid  complications  caused  by  variation 
of  charge,  all  charges  used  in  this  study  were  R5  kg  fired 
at  depths  of  about  100  m.  The  resulting  transfer  function 
is  shown  by  the  dashed  curve  in  Figure  5.  Since  the 
observations  do  not  show  the  high  frequency  ripple  apparent 
in  Figure  5,  we  eliminated  the  ripple  by  applying  a 
minimal  amount  of  absorption.  The  effect  of  absorption 
can  be  treated  by  applying  the  operator  F  (t,  T/Qav) 
defined  by  Carpenter  (1967).  Here  t  is  time,  T  is  the 
total  solid  earth  travel  time  of  the  compressional  wave 
from  source  to  receiver,  and  l/QflV  is  the  average  specific 
dissipation  factor.  The  propagation  operator  used  is 
given  in  Figure  6,  with  T  *  15  sec  and  Qflv  =  500.  Convolving 
this  operator  with  the  transfer  function  produces  the  solid 
curve  in  Figure  5.  This  filtered  transfer  function  was 
used  throughout  the  study. 

Using  the  propagation  operator  in  this  way  is  equivalent 
to  simply  low  pass  filtering  our  synthetic  records.  We  are 
aware  that  while  it  is  operationally  convenient  to  define 
T/Q  as  a  constant,  it  is  an  unrealistic  treatment  of 
absorption  as  it  assumes  the  energy  loss  from  absorption 
does  not  change  with  distance.  The  result  is  that  the 
arrivals  at  short  ranges  are  over-attenuated,  while  those 
at  large  ranges  are  under-attenuated.  Nevertheless,  by 
applying  only  a  small  amount  of  absorption  the  resulting 
amplitude  change  is  not  severe  and  does  not  change  the 
basic  behavior  of  the  amplitude-distance  decay  curves. 

Starting  from  the  travel  time  model  a  perturbed  model 
was  found,  by  trial  and  error,  which  shows  definite 
similarities  in  wave  form  to  the  observed  records;  however, 
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the  refraction  travel  times  now  disagree  slightly.  This 
model ,  F,  is  given  in  Figure  7.  The  densities  and  shear 
velocities  are  assumed  except  for  the  apparent  velocities 
of  3.55  km/sec  and  3.85  km/sec  which  were  obtained  from 
the  travel  time  data.  Before  discussing  the  synthetic 
waveforms  for  this  model  it  is  educational  to  look  at  the 
pressure  response  of  the  model  as  a  function  of  range,  R, 
assuming  a  unit  step  function  pressure  source  (see  Figures 
8  and  9) .  The  pressure  responses  are  the  summation  of  the 
P-mode  generalized  rays  including  multiple  reflections 
within  the  various  layers.  Contributions  of  the  large 
number  of  generalized  rays  that  travel  across  one  layer  or 
more  in  the  S-mode  are  neglected,  although  such  contributions 
are  appreciable  at  short  ranges.  This  subject  will  be 
discussed  in  Part  II.  At  R  ■  14.1  km  the  pressure  response 

is  dominated  by  three  events  which  arrive  at  approximately 

* 

the  same  time.  They  are  responses  from  the  bottom  of 
layers  3,  4,  and  5.  These  three  events  are  quite  strong 
since  all  three  are  near  their  respective  critical  angles, 
hence  the  large  amplitude.  As  the  range  increases  the 
response  from  the  6. 8/8.1  interface,  Pg ,  grows  rapidly.  It 
becomes  the  dominant  event  for  R  >  25  km.  For  larger 
ranges  the  refraction  and  reflection  for  this  interface 
separate  in  time.  As  separation  occurs  the  response  from 
the  bottom  interface,  P^,  grows  supporting  the  refracted 
or  head  wave.  The  peak  following  Pg  is  the  first  internal 
reflection  in  the  oceanic  layer. 

The  synthetic  waveforms  produced  by  convolving  the 
responses  in  Figures  8  and  9  with  the  filtered  transfer 
function  in  Figure  5  are  given  in  Figures  10  and  11.  The 
observations  are  included  for  comparison.  At  the  nearest 
ranges  the  pressure  response  is  essentially  a  step  function 
so  that  when  we  perform  the  convolution  and  take  the 
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derivative  we  reproduce  the  filtered  transfer  function. 

Note  the  similarity  between  Figure  5  and  the  synthetic 
waveforms  at  R  =  14.1  and  18.8  km.  At  larger  ranges  the 
waveforms  of  the  synthetics  are  elongated  by  the  buildup 
of  Pg,  the  Moho  reflection.  The  simplicity  of  the  pressure 
response  at  30.9  and  35.0  km  is  displayed  in  their 
waveforms.  As  we  increase  the  range  the  mantle  head  wave 
begins  to  separate  from  the  Moho  reflection.  This  is  also 
apparent  in  the  observations  especially  at  R  *  47  km. 

However,  the  corresponding  synthetic  record  has  not  reached 
the  same  degree  of  separation.  This  implies  that  the 
critical  range,  that  range  where  the  Moho  reflection  is 
at  critical  angle,  for  the  model  occurs  too  late.  In  fact, 
examining  the  synthetic  with  the  corresponding  observed 
waveforms  we  see  that  the  first  half  of  the  profile.  Figure 
10,  would  agree  much  better  if  the  critical  range  was 
reduced  by  5  km.  This  would  effectively  match  the  synthetic 
waveforms  of  R  «  30.9  and  35.9  km  with  observed  P.  «  27.1 
and  30.9  km.  The  comparison  between  synthetic  and  observed 
waveforms  is  much  better  for  the  second  half  of  the  profile. 
Figure  11,  which  is  controlled  by  the  nature  of  the  crust- 
mantle  transition.  We  will  devote  most  of  this  paper  to 
this  transition.  Adjustments  to  the  upper  crust  where 
converted  shear  waves  are  also  important  will  be  given 
in  Part  II. 

Examining  the  observed  waveforms  of  Figure  11  we  see 
many  variations  in  amplitude  and  shape.  Many  of  these 
variations  are  probably  caused  by  using  a  sinking  explosive 
as  a  signal  generator  with  the  lack  of  exact  positioning 
and  variable  yield.  However,  it  does  appear  that  the 
mantle  headwave  at  separation  is  much  too  large  compared 
to  the  Moho  reflection  to  be  a  true  refraction.  Furthermore, 
it  decays  much  slower  than  the  Moho  reflection.  We  have 
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tried  to  produce  these  effects  by  introducing  the  layer 
in  the  upper  mantle.  We  think  of  this  layer  aa  representing 
a  velocity  gradient  since  replacing  it  by  two  thin  layers, 
maintaining  approximately  the  same  velocity  depth  function, 
produced  about  the  same  synthetic  waveforms.  If  the 
upper  mantle  has  a  low  Q  then  the  gradient  would,  obviously, 
need  to  be  greater.  However,  if  the  gradient  is  increased 
any  appreciable  amount  it  would  be  indicated  on  the  travel 
time  plot  which  it  is  not. 

Two  other  models  were  considered,  one  where  the 
bottom  interface  or  mantle  gradient  was  omitted  and  one 
where  the  mantle  had  a  negative  velocity  gradient;  that 
is  the  bottom  half  space  was  given  the  velocity  7.92  km/sec. 
The  pressure  responses  are  shown  in  Figure  12  and  the 
synthetic  waveforms  in  Figure  13  for  R  =  84  km.  We  include 
the  model  containing  t'\e  velocity  increase  for  comparison. 
Figure  13  suggests  that  if  the  upper  mantle  in  this 
region  has  a  constant  velocity  it  would  be  very  difficult 
to  identify  a  head  wave  at  ranges  greater  than  about  80  km 
assuming  the  same  source  and  receiver.  The  case  of  a 
negative  gradient  would  damp  out  the  head  wave  even  faster. 

Finally  we  consider  the  sharpness  of  the  Moho 
transition.  Examining  Figure  11  we  find  a  relatively  quiet 
zone  between  the  head  wave  and  Moho  reflection.  If  there 
is  a  significant  transition  layer  between  the  crust  and 
the  mantle  its  arrivals  would  manifest  themselves  in  this 
zone.  To  illustrate  this  point  we  constructed  a  model 
containing  a  one-kilometer  transition  layer  with  velocity 
of  7.5  km/sec.  The  results  at  R  »  84  are  given  in  Figure 
14.  Since  the  comparison  with  observations  is  obviously 
poor  we  conclude  that  the  transition  is  sharp,  probably 
occurring  over  much  less  than  a  kilometer  of  thickness  at 
this  particular  station. 


-87- 


COMPARISON  OF  THEORETICAL  AND  OBSERVED  AMPLITUDES 

The  amplitude-distance  decay  curve  for  the  synthetic 
records  was  measured  using  the  same  method  as  for  the 
observed  records.  The  comparison  is  given  in  Figure  15. 
While  the  comparison  is  not  perfect,  the  decay  curve  for 
the  synthetic  first  arrivals  has  the  same  general 
characteristics  as  that  of  the  observed  records.  However, 
the  observed  amplitudes  do  decay  more  rapidly  with  distance 
than  do  the  synthetics  over  the  first  50  km.  This  is 
probably  caused  by  not  having  the  correct  upper  crustal 
model.  After  studying  the  ratio  of  the  primary  waves  to 
converted  shear  wave.  Part  II,  we  found  it  necessary  to 
reduce  the  P-velocities  of  the  upper  crust.  This 
adjustment  produced  a  faster  amplitude  decay  with  range. 
Nevertheless,  the  large  observed  amplitude  at  R  ■  18.8, 
as  well  as  some  of  the  rather  erratic  amplitude  behavior 
at  large  ranges,  is  not  explained.  It  does  not  appear 
that  a  plane,  homogeneous,  isotropic  layered  model  can 
explain  such  jumps. 

CONCLUSION 

We  have  shown  that  the  velocity  must  increase  rapidly 
at  the  sea  bottom  to  produce  the  large  amplitudes  of  the 
first  arrivals  at  small  range.  At  larger  ranges  the 
existence  of  Moho  reflections  indicates  the  sharpness  of 
the  crust-mantle  transition.  In  this  region  the  transition 
zone  is  probably  much  less  than  a  kilometer  in  thickness. 

The  slow  decay  of  mantle  head  waves  with  range  is  explained 
by  a  positive  velocity  gradient  in  the  upper  few  kilometers 
of  the  mantle.  To  what  depth  this  gradient  is  maintained 
cannot  be  answered  by  this  study  due  to  limited  observations. 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 

Figure  3 


Figure  4 
Figure  5 

Figure  6 
Figure  7 
Figure  8 


Sample  recording  showing  prominent  arrivals. 

The  upper  three  traces  are  high  gain,  low 
frequency  (3-20Hz)  outputs  from  three  different 
hydrophones.  The  center  trace  is  of  inter¬ 
mediate  frequency  (20-200Hz).  The  lower  three 
traces  are  low  frequency  outputs  attenuated 
in  magnitude  by  20  db  compared  to  the  upper 
three. 

Reduced  travel  time  plot  and  layered  model 
interpretation  for  refraction  profile  SH31. 

Amplitude  decay  with  distance  for  profile  SH31. 
Second  arrivals  correspond  to  those  arrivals 
in  Figure  2  which  lie  on  or  near  the  6.94  km/sec 
line  and  are  at  distances  greater  than  50  km. 
These  are  a  combination  of  head  waves  from  the 
oceanic  layer  and  reflections  from  the  Moho- 
rovicic  discontinuity. 

Amplitude  decay  with  distance  for  profile  SH11. 

Transfer  function  before  and  after  low-pass 
filtering. 

Propagation  Operator. 

Oceanic  crustal  model  F. 

Pressure  response  versus  time  after  onset 
assuming  a  step  function  pressure  source.  P^ 
indicates  the  reflection  time  from  the  bottom 
layer  i. 
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Figure  9  Continuation  of  Figure  8.  Pg  indicates  the 
Moho  reflection  time. 

Figure  10  Synthetic  waveforms  based  on  the  pressure 
responses  of  Figure  8. 

Figure  11  Synthetic  waveforms  based  on  the  pressure 
responses  of  Figure  9. 

Figure  12  Pressure  responses  for  a  crustal  model  over- 
lying  a  mantle  with  positive  gradient, 
negative  gradient  and  constant  velocity. 

Figure  13  Synthetic  waveforms  based  on  pressure 
responses  of  Figure  12. 

Figure  14  Pressure  response  and  synthetic  waveform  for 
a  modal  containing  a  transition  layer  at 
range  84  km. 

Figure  15  Comparison  of  relative  amplitude  decay  with 

distance  for  observed  records  from  profile  SH31 
and  for  theoretically  calculated  synthetic 
records  for  model  F. 
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Ill .2-2:  Converted  Shear  Waves 


In  a  recent  paper  by  the  authors  part  1,  (Helmberger , 
Morris,  1968)  synthetic  primary  waves  for  some  oceanic 
crustal  models  were  compared  with  a  profile  of  observations. 
At  short  ranges,  R  <  50  km,  a  rather  large  second  arrival 
is  observed  which  has  the  apparent  velocity  of  a  shear 
wave,  see  Figures  (2)  and  (3).  We  will  label  such  arrivals 
converted  shear  waves,  waves  that  travel  through  the  ocean 
as  primary  waves  and  in  the  solid  crust  as  shear  waves. 

It  is  found  that  these  generalized  rays  have  a  behavior 
similar  to  primary  waves.  That  is  the  summation  of  rays 
leads  to  a  large  arrival  which  increases  its  velocity  as 
a  function  of  range. 

The  model  studied  in  part  1  agreed  well  at  large 
ranges  (50  to  112  km)  where  the  waveforms  are  dominated 
by  the  structure  near  the  crust-mantle  transition.  The 
model  indicated  a  sharp  crust-mantle  boundary  with  a  small 
positive  velocity  gradient  in  the  upper  mantle.  However, 
this  model  did  not  agree  very  well  at  short  ranges  where 
the  critical  Moho  reflection  of  the  model  occurs  at  about 
5  kms  greater  range  than  the  observed.  Before  computing 
converted  shear  waves  we  will  adjust  the  crustal  structure 
to  rectify  this  apparent  disagreement. 
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Mode  Conversion 

One  oi  tl.e  disadvantages  of  using  generalized  ray  the¬ 
ory  in  syntl  etic  seisuogram  generation  is  the  troublesome 
problem  of  conversion  of  wa  »e  types  at  interfaces.  This 
would  not  be  such  a  difficulty  if  one  was  not  forced  to 
use  i^any  layers  to  approximate  crustal  structures.  However, 
it  appears  that  many  generalized  rays  containing  mixed  code 
paths  can  be  neglected  at  certain  ranges.  In  a  recent  paper 
(Hel  iher^er ,  1968)  a  procedure  was  presented  for  calculating 
the  response  for  a  sandwiched  Layer,  including  nixed  node 
paths.  The  numerical  results  for  a  one  layer  crustal  model 
is  given  in  Fig.  1  where  the  transfer  function  of  part  1 
was  used  as  it  will  be  throughout  this  study.  The  seismo¬ 
grams  at  R-19  and  32  km  contain  the  first  and  second  set  of 
generalized  rays.  Only  the  first  set  is  included  at  R-25  km. 
It  is  clear  from  Fig.  1  that  the  first  set  yields  the  prin¬ 
cipal  contributions  at  the  larger  ranges.  The  results  show 
that  P  to  S  mode  conversion  is  the  strongest  at  shorter 
ranges  and  drops  off  rapidly  approaching  the  P-mode  criti¬ 
cal  angle.  Fig.  1  should  be  compared  with  Fig.  2  which  dis¬ 
plays  part  of  the  profile  of  observations  studied  in  part  1. 
Fig.  3  is  along  the  same  profile  but  where  two  pound  charges 
were  used  to  generate  the  seismic  signals.  Note  the  differ¬ 
ence  in  the  principal  period  of  the  waveforus  of  Fig.  3  com¬ 
pared  to  those  of  Fig.  2.  This  is  caused  by  the  variation 
in  bubble  size  which  is  predictable.  The  comparison  of  Fig. 

1  with  Fig.  2  and  3  suggest  that  PP  and  SS  are  idealized 
and  Gj  respectively  at  the  larger  ranges.  This  identifica¬ 
tion  will  be  supported  by  this  study.  Ond  of  the  largest 
discrepancies  is  the  ratio  of  PP  to  SS  as  compared  to  G^ 
and  G2»  This  is  caused  by  the  large  P-velocity  contrast  at 
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tl  e  ocean  bottoc.  interface.  By  reducing  the  upper  crustal 
velocity  the  ratio  of  PP  to  SS  will  be  enhanced,  also  the 
a  plitude  of  Pr ,  now  a  diffracted  wave  caused  by  the  gra¬ 
dient,  will  be  greatly  increased  ac  the  nearest  ranges. 

The  travel  tivfie  between  Pr  and  PP  will  also  be  reduced  and, 
as  we  will  see,  ttiese  two  events  will  interfere  approaching 
the  PP  critical  angle.  Note  the  disturbances  in  and  G9 
at  about  R“25  kr...  Examining  the  top  three  recordings  in 
Fig.  3  it  appears  that  the  signals  between  G^  and  G could 
well  le  rode  conversions  since  they  decay  rapidly  with  range. 
It  would  Le  gratifying  to  show  this  f on  ally  by  sur.iv.iing 
all  possible  generalized  rays  but  t  iis  would  be  very  tie c 
consuiin0.  We  will  neglect  these  ode  conversions  in  this 
study  and  compute  G^  based  oh  P-  iode  and  G ^  based  on  S-uode 
exclusively.  That  is  we  will  allow  changes  only  in  i.:ode 
crossing  the  ocean  bottom,  interlace. 
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Crustal  Model 

The  structure  near  the  crust -uant le  transition  was 
adequately  determined  in  part  1  which  we  will  incorporate 
into  our  new  model.  Suppose  we  reinterpret  the  travel  time 
plot  given  in  Fig.  4  by  adding  layers  in  the  upper  crust. 

The  added  criterion  being  that  the  new  model  has  its  cri¬ 
tical  Mono  reflection  at  about  ?.<}  lo  •  and  that  the  Molio 
reflect  Lon ,  at  large  ranges,  1  as  the  proper  lag  behind  the 
lead  wave.  The  model  candidates  can  be  easily  checked 
aGaiust  t.  c  above  criterion  by  applying  the  first-motion 
approximations  discussed  by  Helaberger  (1968).  A  model 
passing  these  tests  is  given  in  Fig.  4,  found  by  trial  and 
error.  The  model  contains  the  same  gradient  in  the  upper 
mantle  as  discussed  in  Part  1,  that  is  .05  km/sec  per  kn. 
Synthetic  seismograms  constructed  exclusively  from  P-mode 
paths  are  shown  in  Fig.  5  which  should  be  compared  with 
Figures  2  and  3.  At  the  shortest  ranges  the  summation  of 
reflections  from  the  upper  crustal  interfaces  produce  the 
principal  motion  which,  may  be  interpreted  as  the  direct  wave 
turned  back  by  the  velocity  gradient.  Interference  is  pro¬ 
duced  when  this  event  approaches  the  Moho  reflection  at 
larger  ranges  as  dramatically  indicated  at  R-22.7  in  Figures 
2  and  5.  The  same  phenomeno  occurs  in  G£  at  R“27.1.  As 
we  increase  the  range  the  synthetic  seismograms  approach 
those  generated  in  part  1  since  the  lower  structure  is  simi¬ 
lar. 
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Converted  Shear  Waves 

The  S-velocity  model  was  determined  by  applying  the 
standard  refraction  interpretation  technique  but  where 
tlie  layer  thicknesses  are  already  assigned.  The  resulting 
model  is  not  as  well  determined  as  the  P-velocity  model  since 
the  observations  are  limited  to  40  km.  The  model  used  in 
this  study  is  given  in  Fig.  4.  A  profile  of  pressure  responses 
based  on  this  model  is  displayed  in  Fig.  6,  assuming  a  step 
function  source.  The  responses  are  the  summation  of  genera¬ 
lized  rays  including  multiple  reflected  paths  within  the 
various  layers.  The  long  period  precursor  is  the  Rayleigh 
wave  which,  of  course,  is  not  excited  by  our  source  function. 
The  crustaL  arrivals  add  up  to  produce  the  large  negative 
peak  Sg  at  R-14.1  km  which  can  be  interpreted  as  the  direct- 
wave  turned  around  by  the  velocity  -radical.  As  we  go  to 
larger  ranges  the  Moho  reflection  builds  up  denoted  by  S^. 

These  two  arrivals  interact  as  they  approach  each  oti  er. 
Convolving  the  pressure  response  wici  the  derivative  of  the 
transfer  function  with  respect  to  tii  c  produces  the  synthe¬ 
tic  seismograms  plotted  in  Fig.  b.  Comparing  Fig.  5  with 
the  corresponding  observations  displayed  in  Fig.  2  we  find 
a  reasonably  good  fit  except  that  the  observed  shear  waves 
appear  to  die  out  rather  abruptly.  We  car.  not  explain  this 
feature  present ly.  Furikoru.ore ,  liic  P-  ode  generalized  ray 
that  is  once  reflected  in  the  ocean  layer  arrives  on  top  of 
the  shear  wave  at  about  42  km  which  complicates  the  situation. 

A  number  of  other  stations  are  row  be  in  ,  studied  to  help  cla¬ 
rify  this  discrepancy  and,  if  possible,  identify  the  ant le 
shear  head  wa \re . 


-112- 


| 

Figure  Captions 

Figure  1 

Figure  ? 

Figure  3 

Figure  4 
Figure  5 

Figure  6 


One  layer  oceanic  crustal  model  and  syn- 
thetic  responses .  Velocities  in  ta/sec, 
thicknesses  in  km. 

Observed  recordings  frou  a  protiLe  of 
d5  kg  charges. 

Observed  recordings  frou  a  profile  of 
2  pd  charges.  AH  recordings  on  the  same 

amplitude  scale. 

Reduced  travel  tine  plot  and  layered  .odel 

interpretation . 

Synthetic  seismograms  from  a  Pr°fij-e  °f 
83  kg  charges  assuming  the  model  given 

Fig.  4. 

Pressure  response  versus  time  after  onset 
assuming  a  step  function  pressure  sou 
Sq  is  the  summation  of  upper  layer  responses 
„„.j  Ci*  Hie  Moho  response. 
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The  research  efforts  of  the  research  associates  supported  by  this 
program  are  described.  These  studies  are  divided  into  two  groups: 

I.  Continental  size  arrays  and  their  application,  and  II.  Theoretical 
seismograms  and  earth  structure. 

A  continental  size  array  can  be  realized  by  considering  individual 
seismic  stations  as  array  elements.  A  section  of  North  America,  where 
LRSM  stations  exist  in  sufficient  density,  is  considered  as  a  very 
large  aperture  array  to  supplement  LASA.  Travel  times  and  dT/dA 
values  for  core  phases  (PKP,  PKIKP)  are  measured  with  this  array  and 
the  results  are  interpreted.  A  technique  for  direct  measurement  of 
dT/dA  is  also  described. 

The  theoretical  seismograms  for  P  and  S  waves  are  computed  for  a 
model  of  point  source  in  layered  medium  and  matched  with  observed 
records.  With  this  technique,  the  structures  are  determined  for 
oceanic  crust  using  marine  refraction  profiles,  and  for  the  upper 
mantle  utilizing  seismograms  from  NTS  nuclear  explosions. 


J473 


Unclassified 


Srt’iiritv  Cl.issiln  «iti.»n 


Unclassified 


Seismology 
Earth’s  core 
Arrays 

Synthetic  seismograms 
Mantle  structure 
Travel  times 


HOLE  *  r 


Unclassified 


Security  Classification 


